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We review tiie fictitious integrable system approacli which predicts dynamical tunneling rates 
from regular states to the chaotic region in systems with a mixed phase space. It is based on the 
introduction of a fictitious integrable system that resembles the regular dynamics within the regular 
island. We focus on the direct regular-to-chaotic tunneling process which dominates, if nonlinear 
resonances within the regular island are not relevant. For quantum maps, billiard systems, and 
optical microcavities we find excellent agreement with numerical rates for all regular states. 



o 

(N 

> 
O 

:z; 

(N 



Q 
U 



> 
00 

O 

a\ 
o 
o 



PACS numbers: 05.45.Mt, 03.65.Sq, 03.65.Xp 



I. INTRODUCTION 

Tunneling of a particle is one of the central manifes- 
tations of quantum mechanics. The prototypical exam- 
ple is the tunneling escape from a one-dimensional po- 
tential well through an energy barrier. While classically 
the particle is confined for all times, quantum mechani- 
cally the probability inside the well decays exponentially, 
exp(— 7t), where 7 is the tunneling rate. It depends on 
the width and the height of the barrier and can be pre- 
dicted, e.g., using WKB theory [jy]. Tunneling vanishes 
in the semiclassical limit, where typical classical actions 
are large compared to Planck's constant. 

Tunneling not only occurs for potential barriers but 
whenever the corresponding classical system consists of 
dynamically disconnected regions in phase space, which 
has been termed dynamical tunneling p]. It occurs in 
Hamiltonian systems which typically have a mixed phase 
space. Here regions of regular motion, the so-called reg- 
ular islands, and regions of chaotic motion, the so-called 
chaotic sea, coexist. While the classical motion is con- 
fined to any of these regions, quantum mechanically they 
are coupled by dynamical tunneling. In particular the 
fundamental process of regular-to- chaotic tunneling de- 
scribes the exponential decay of a wave packet initially 
localized in the regular island to the chaotic sea. The 
same coupling also leads to tunneling from the chaotic 
sea to the regular island. 

Dynamical tunneling also affects the structure of eigen- 
states of systems with a mixed phase space. According 
to the semiclassical eigenfunction hypothesis the 
eigenstates are concentrated either in the regular islands 
or in the chaotic sea. Away from the semiclassical limit 
this classification still holds approximately such that the 
corresponding eigenstates are called regular or chaotic. 
However, each eigenstate has contributions in the other 
regions of phase space, due to dynamical tunneling. 

Dynamical tunneling in a mixed phase space has 
been studied theoretically [6|436l and experimentally, 
e.g. , in cold atom systems j37h39j , microwave billiards 
[33 . liol Idlj , and semiconductor nanostructures [i^l ■ It is 
of current interest for, e.g., eigenstates affected by flood- 



ing of regular islan ds |4 3'-^i6l , emission properties of op- 
tical microcavities [SJ, [47, - 49.] , and spectral statistics in 
systems with a mixed phase space [50|-|53|. 

The concept of chaos-assisted tunneling was intro- 
duced in Refs. (sl-fioj. It occurs, e.g., in systems with two 
symmetry-related regular islands surrounded by chaotic 
motion in phase space. There it was observed that the 
energy splittings /S.E between two symmetry related reg- 
ular states are typically drastically enhanced due to the 
appearance of chaotic states, compared to an integrable 
situation. Chaos- assisted tunneling consists of two pro- 
cesses: a regular-to-chaotic tunneling step from a regu- 
lar torus of one island to the chaotic sea and a chaotic- 
to-regular tunneling step from the chaotic sea to the 
symmetry-related torus. The energy splittings IS.E show 
strong fluctuations [sl-fiol. Is^ ] under variation of exter- 
nal parameters, as the distance between the energies 
of the regular doublet and the close-by chaotic states 
varies. This was also observed for optical microcavi- 
ties and microwave billiards [13, Ej . In contrast to 
the energy splittings IS.E the regular-to-chaotic tunneling 
rates 7 describe the average coupling to the chaotic sea, 
7 — (AE) /h, as shown in Sec. Ill Dl and therefore do not 
fluctuate. 

In the regime, h < ^rog, in which the Planck constant 
h is smaller but of the same order as the area Aj-cg of 
the regular island, quantum mechanics is not affected by 
fine-scale structures in phase space, such as nonlinear res- 
onances or the hierarchical transition region. Hence, the 
direct regular-to- chaotic tunneling process dominates. In 
Ref. a qualitative argument for the tunneling rates to 
behave exponentially as 7 cx exp{—B/h) was given. One 
may write B = CA^-cg, where the non- universal constant 
C has been calculated for various situations in the semi- 
classical limit: For a weakly chaotic system C = 2tt [sO] 
was found and C > 2tt for rough nanowires [H, [s^ . The 
approach in Ref. [2^ leads semiclassically to C = 2 — In 4 
which was corrected to C = 3 — In 4 [57j. However, this 
prediction does not describe generically shaped regular 
islands as it uses a transformation of the regular island 
to a harmonic oscillator. 

In order to find a quantitative prediction of direct 
regular-to-chaotic tunneling rates for generic islands 
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we introduced the fictitious integrable system approach 
[3^ [ssl . It rehes on the decomposition of Hilbert space 
into a regular and a chaotic subspace by means of a fic- 
titious integrable system [2^ IsO] ■ We require that its 
dynamics resembles the regular motion in the originally 
mixed system as closely as possible and extends it beyond 
the regular region. This leads to a tunneling formula in- 
volving properties of this integrable system as well as 
its difference to the mixed system under consideration. 
It allows for the prediction of tunneling rates from any 
quantized torus within the regular island. This approach 
was applied to quantum rnaps [32j . billiard systems [ssj , 
and optical microcavities |34| . 

In the semiclassical regime, h ^ ^rcg, the fine-scale 
structures of the classical phase space can be resolved by 
quantum mechanics. In particular, nonlinear resonances 
cause an enhancement of the regular-to-chaotic tunneling 
rates, which has been termed resonance-assisted tunnel- 
ing [l8l - l2l| . It leads to characteristic peak and plateau 
structures in the tunneling rates as observed for near in- 
tegrable systems [3, [l^ , mixed quantum maps [1^ |2l| , 
periodically driven systems 22, 23], quantum accelerator 
modes [3d| . and for multidimensional molecular systems 
(26I [27} . Quantitatively, however, deviations of several 
orders of magnitude to numerical rates appear, especially 
in the experimentally accessible regime where nonlinear 
resonances become relevant for tunneling. Recently, it 
was shown that a combination of the direct regular-to- 
chaotic tunneling mechanism and the resonance-assisted 
tunneling mechanism leads to a theory which quantita- 
tively predicts tunne ling rates from the quantum to the 
semiclassical regime [35| . For the application of this uni- 
fied theory it is essential to determine the direct regular- 
to-chaotic tunneling rates. 

In this paper we review the fictitious integrable system 
approach for the prediction of direct regular-to-chaotic 
tunneling rates. The approach is derived in Sec. |lll 
Numerical methods for the determination of tunneling 
rates are presented in Sec. IIIII It is then shown how 
the approach can be applied analytically, semiclassically, 
and numerically to quantum maps in Sec. IIVI and two- 
dimensional billiard systems in Sec. [V] 



II. DYNAMICAL TUNNELING AND THE 
FICTITIOUS INTEGRABLE SYSTEM 
APPROACH 

We consider systems with a mixed phase space, in 
particular two-dimensional quantum maps and billiards. 
Their phase space is divided into regions of regular dy- 
namics and regions of chaotic dynamics. We focus on 
the fundamental situation of just one regular island em- 
bedded in the chaotic sea. Fig. [IJa). At the center of 
the island one has an elliptic fixed point, which is sur- 
rounded by invariant regular tori. For such systems the 
semiclassical eigenfunction hypothesis implies that 
in the semiclassical limit the eigenstates can be classified 



as either regular or chaotic, according to the phase-space 
region on which they concentrate. In order to understand 
the behavior of eigenstates away from the semiclassical 
limit one has to compare the size of phase-space struc- 
tures with Planck's constant h. We discuss this exem- 
plarily for quantum maps (soj which are described by a 
unitary time-evolution operator U on a Hilbert space of 
finite dimension N . Here we introduce the semiclassical 
parameter hes = h/A = 1/N as the ratio of Planck's 
constant h to the area A of phase space. The eigenstates 
\ipn) and quasi-energies ipn of U are determined by 



(1) 



The so-called regular states are predominantly concen- 
trated on tori within the regular island and fulfill the 
Bohr-Sommerfeld type quantization condition 



pdq^ hcs m + 
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0, 



1. (2) 



For a given value of hes there exist iVreg of such regular 
states, where TV^g = [^reg/^off+l/2j and A^g from now 
on is the dimensionless area of the regular island. The 
chaotic states mainly extend over the chaotic sea. Note, 
that for systems with a large density of chaotic states the 
regular states may disappear and chaotic states fiood the 
regular island 0, HI] . 

An important consequence of a finite hes in systems 
with a mixed phase space is dynamical tunneling. It 
couples the regular island and the chaotic sea, which 
are classically separated. Hence, the regular and chaotic 
eigenfunctions of U always have a small component in the 
other region of phase space, respectively, see Fig.[TIb,c). 

The coupling of the regular and the chaotic phase- 
space regions can be quantified by tunneling rates 7^ 
which describe the decay from the mth regular torus to 
the chaotic sea. To define these tunneling rates one can 
consider a wave packet started on this mth quantized 
torus in the regular island which is coupled to a contin- 
uum of chaotic states, as in the case of an infinite chaotic 
sea. Its decay e"'*'"' is characterized by a tunneling rate 
7m- For systems with a finite phase space this expo- 
nential decay occurs at most up to the Heisenberg time 
th = hgff/ Ach, where A^h is the mean level spacing of the 
chaotic states. Alternatively, the tunneling rates can be 
obtained from lifetimes of resonances in a corresponding 
open system, e.g., by adding an absorbing region some- 
where in the chaotic component, see Sec. IIII Al 

In the regime, heS ^ ^reg, where hcs is smaller 
but comparable to the area ^rog of the regular island, 
the rates 7™ are dominated by the direct regular-to- 
chaotic tunneling mechanism, while contributions from 
resonance-assisted tunneling are negligible. We concen- 
trate on systems where additional phase-space structures 
within the chaotic sea are not relevant for tunneling. In 
the following we derive a prediction for the direct regular- 
to-chaotic tunneling rates using the fictitious integrable 
system approach ^32| . 
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FIG. 1. (Color online) (a) Illustration of the mixed classical 
phase space corresponding to a quantum map U together with 
the Husimi representation of (b) a regular and (c) a chaotic 
eigenstate of U which both have a small component in the 
other region, (d) Illustration of the classical phase space of 
the fictitious integrable system l7rog. (e) Eigenstates |i/'rcg) 
of f/reg are purely regular, while (f) the purely chaotic states 
\ipch) extend in the chaotic region of phase space. 



A. Derivation 

In order to find a prediction for the direct regular-to- 
chaotic tunneling rates we decompose the Hilbert space 
of the quantum map U into two parts which correspond 
to the regular and chaotic regions. While classically 
such a decomposition is unique, quantum mechanically 
this is not the case due to the uncertainty principle. 
We find a decomposition by introducing a fictitious in- 
tegrable system U^eg- Related ideas were presented in 
Refs. [^,[2^ [3^. The fictitious integrable system has to 
be chosen such that its dynamics resembles the classical 
motion corresponding to U within the regular island as 
closely as possible and continues this regular dynamics 
beyond the regular island of U, see Fig. [Hd). The eigen- 
states iV'reg) of Ureg, C/reglV'reg) = e"^"HV'reg)> are purely 
regular in the sense that they are localized on the mth 
quantized torus of the regular region and continue to de- 
cay beyond this regular region, see Fig. [TJe) . This is the 
decisive property of IV'reg) which have no chaotic admix- 
ture, in contrast to the predominantly regular eigenstates 
of U, see Fig. [IJb). The explicit construction of Urcg is 
discussed in Sec. IIIBl 

With the eigenstates iV'reg) ^reg we define a projec- 
tion operator 

Preg:= E l«(^"gl, (3) 

using the first iVreg eigenstates of C/reg which approxi- 
mately projects onto the regular island corresponding to 



U. The orthogonal projector 

Pch := 1 - Prog (4) 

approximately projects onto the chaotic phase-space re- 
gion. These projectors, Prcg and Pch, define our decom- 
position of the Hilbert space into a regular and a chaotic 
subspace. 

Introducing a basis |'0ch) in the chaotic subspace we 
can write Pch = J2ch lV'ch)(V'ch|- Here we sum over all 
^ch = N — Ni-cg states \ipch), which we call purely chaotic 
states, see Fig.[TJf) for an illustration. The coupling ma- 
trix element Vch,m between a purely regular state IV'rog) 
and any purely chaotic state IV'ch) is 

«ch,™ = (^ch|C/|^"g>- (5) 

From this the tunneling rate is obtained using a dimen- 
sionlcss version of Fermi's golden rule, see Appendix Rl 

7m = ^ bch.mP, (6) 
ch 

where the sum is over all chaotic basis states \tpch) and 
thus averages the modulus squared of the fluctuating ma- 
trix elements Vch,m- Here we apply Fermi's golden rule 
in the case of a discrete spectrum, which is possible if 
one considers the decay e"'*'"* up to the Heisenberg time 
TH = Ks/ Ach only. 

Inserting Eq. ([5]) into Eq. ([5]) we obtain 

7m = \\Pci.umg)r = 11(1 - Prog)c/i^rcg)f (7) 

as the basis of all our following investigations. It allows 
for the prediction of tunneling rates from a regular state 
localized on the mth quantized torus to the chaotic sea. 
Equation ([7]) agrees with the intuition that the tunnel- 
ing rates are determined by the amount of probability 
that is transferred to the chaotic region after one appli- 
cation of the time evolution operator U on IV'reg)- We 
want to emphasize that Eq. ([7]) essentially relies on the 
chosen decomposition of Hilbert space determined by the 
fictitious integrable system J/rog- A similar expression for 
the tunneling rates was obtained from a phenomenolog- 
ical Hamiltonian in Ref. [s^. Note, that the tunneling 
rate for the inverse process of tunneling from the chaotic 
sea to the mth regular torus is also given by Eq. ([7]) but 
with an additional prefactor of 1 /Nch due to the different 
density of final states [46j . 

1. Approximation for very good Ureg 

In cases where one finds a fictitious integrable system 
C/rcg which resembles the dynamics within the regular 
island of U with very high accuracy, Eq. (O can be ap- 
proximated as 

7m«ll(t^-C/reg)|V'"g)lP, (8) 

using PrcgC/|?/'r"og) ~ -Prcg C/rog IV'reg) = C^rog|Vr"g>- Instead 
of the projector Pch in Eq. ([7]) the difference U —U^cg en- 
ters in Eq. ([5]) . This allows for a semiclassical evaluation, 
which is presented in Sec. Ill CI 
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2. Approximation for non-orthogonal chaotic states 

If one constructs chaotic states \ipch) from random 
wave models Q , they will not be orthogonal to the purely 
regular states IV'reg)- this case we construct orthonor- 
malized states 

IV'ch) :=c(ll-Preg)|V^ch) (9) 

with normalization c. We find for the coupling matrix 
elements, Eq. JS]), 

«(^eh|C/-C/rcg|V'"g>, (11) 

where we use the approximations c w 1 and again 
Preg?7|V'™g) ~ C^rcg|V'"g)- Equation ([11]) can now be in- 
serted into Eq. ([6]), leading to 

7m«ll/^ch(t/-C/rcg)|^rcg)ll' (12) 
with Pch = I]ch IV'ch)(^ch|- 

3. Application to billiards 

Two-dimensional billiard systems, which we consider 
in Sec. |Vl are given by the motion of a free particle 
of mass M in a domain il with elastic reflections at 
its boundary dQ. Quantum mechanically they are de- 
scribed by a Hamilton operator H. The fictitious inte- 
grable system approach can also be applied to billiards: 
We use a fictitious integrable system i/rcg and its eigen- 
states ■(/'™g (q), characterized by the two quantum num- 
bers m and n. We start from a random wave model Q 
for the chaotic states V'ch(9) which are not orthogonal 
to the purely regular states. Using the approximation 
for non-orthogonal chaotic states we obtain in analogy 
to Eq. (HID 

Kh,m„ = J d\Mq){H - i7,eg)^™(q) (13) 
n 

for the coupling matrix element between a purely regular 
state with quantum numbers (m, n) and different chaotic 

states 'ipchiq)- The tunneling rate Tmn is obtained using 
Fermi's golden rule, Eq. (IA1[) . 

Tmn = — (|Vch.mn| ) Pch ~ , , , (|K:h.mK| ) (14) 

n 4M 

where we average over the modulus squared of coupling 
matrix elements Vch.rnn between one particular purely 
regular state and different chaotic states of similar en- 
ergy. The chaotic density of states pch is approximated 
by the leading Weyl term p^h ~ Achfi^/ (SttM) in which 
Ach denotes the area of the billiard times the chaotic 
fraction of phase space. This expression for p^h follows, 
e.g., from counting the number of Planck cells in the 
chaotic part of phase space 



B. Fictitious integrable system 

The most difficult step in the application of Eqs. ([7]) 
and ()14|) to a given system is the determination of the 
fictitious integrable system H^cg- Its dynamics should 
resemble the classical motion of the considered mixed 
system within the regular island as closely as possible. 
As a result the contour lines of the corresponding in- 
tegrable Hamiltonian i/icg, Fig. [TJd), approximate the 
KAM-curves of the classically mixed system. Fig. [Ija) , in 
phase space. This resemblance is not possible with arbi- 
trary precision as the integrable approximation for exam- 
ple does not contain nonlinear resonance chains and small 
embedded chaotic regions. Moreover, it cannot account 
for the hierarchical regular-to-chaotic transition region 
at the border of the regular island. Similar problems ap- 
pear for the analytic continuation of a regular torus into 
complex space due to the existence of natural boundaries 
[^ [lll - [Ti . [l8l - [20i . However, for not too small hes, where 
these small structures are not yet resolved quantum me- 
chanically, an integrable approximation with finite ac- 
curacy turns out to be sufficient for a prediction of the 
tunneling rates. 

In addition the integrable dynamics of i?rcg should ex- 
trapolate smoothly beyond the regular island of H. This 
is essential for the quantum eigenstates of i?reg to have 
correctly decaying tunneling tails. According to Eq. ([7]) 
they are relevant for the determination of the tunneling 
rates. While typically tunneling from the regular island 
occurs to regions within the chaotic sea close to the bor- 
der of the regular island, there exist other cases, where it 
occurs to regions deeper inside the chaotic studied 
in Ref. [30]. Here i?rcg has to be constructed such that 
its eigenstates have the correct tunneling tails up to this 
region, see Sec. IIV A3I 

For quantum maps we determine the fictitious inte- 
grable system in the following way: We employ classical 
methods, see below, to obtain a one-dimensional time- 
independent Hamiltonian H,-cg{q,p) which is integrable 
by definition and resembles the classically regular motion 
of the mixed system. After its quantization we obtain 
the regular quantum map U^cg = e"'^'"*/''''" with corre- 
sponding eigenfunctions IV'reg)- F'^^ the numerical evalua- 
tion of Eq. (0) we use Peh = 1-Prog = 1-E IV'"g)('/'"gl, 
where the sum extends over m = 0, 1, . . . , N^cg — 1. 

Now we discuss two examples for the explicit construc- 
tion of H^cg- Note, that also other methods, e.g., based 
on the normal- form analv sis (63 , |6]| or on the Campbell- 
Baker-Hausdorff formula [62| can be employed in order 
to find Hj-cg. For the example systems considered in this 
paper, however, they show less good agreement. 

1. Lie-transformation method 

One approach for the determination of the ficti- 
tious integrable system for quantum maps is the Lie- 
transformation method 63'|. It determines a classical 
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Hamilton function 

K 

H^,^{q.p) ^Y.^'hi{q,p) (15) 

1=0 

as a power series in the period of the driving r, see 
Fig. [2ja) and Ref. jl^ for examples. Typically, the order 
of the expansion K can be increased up to 20 within rea- 
sonable numerical effort. The Lie-transformation method 
provides a regular approximation H^-cg which interpo- 
lates the dynamics inside the regular region and gives 
a smooth continuation into the chaotic sea. At some 
order K the series typically diverges due to the nonlin- 
ear resonances inside the regular island. For strongly 
driven systems, such as the standard map at k > 2.5, the 
Lie-transformation method is not able to reproduce the 
regular dynamics of t/, see Fig. [Hb). 

2. Method using the frequency map analysis 



for each torus in the regular region of phase space [58|. 
This formula follows from Hamilton's equations of motion 
and A{u) = § p{q, u) dq. Finally, we find the fictitious in- 
tegrable system Hrcg by a two-dimensional extrapolation 
of the energies to the whole phase space with 

K 

H^.g{<l,p)= E /ife,/e'""='e2-'f (17) 

k,l=-K 

using periodic basis functions up to the maximal order K. 
An example is shown in Fig.[3l^a). Note, that the result- 
ing -ffrog shows a reasonable behavior beyond the regular 
island of H only for small values of K. For too large 
orders K one observes that 7?reg oscillates in this region, 
see Fig. [3)Jb). This would lead to purely regular states 
IV'reg) with incorrect tunneling tails beyond the regular 
island of J7, resulting in wrong predictions of tunneling 
rates with Eq. ^ . 



An alternative method is applicable even to strongly 
driven one-dimensional systems. In order to determine 
Hrog{q,p) we associate to each torus within the regu- 
lar region of U an energy. This information for in- 
dividual tori is then extrapolated to the entire phase 
space. To this end we consider a straight line {q{u) , p{u)) , 
parametrized by u, from the center u = of the regular 
island to its border with the chaotic sea. Each torus of 
the map crosses this line at some value u and using the 
frequency map analysis [6^ we compute the enclosed area 
A(u) and the rotation number lyiu). Using a polynomial 
interpolation of these functions we calculate an energy 

u 

E{u) = I du'i.(u')^4^ (16) 
J du' 






FIG. 2. (Color online) Application of the Lie-transformation 
method, (a) Orbits (thick gray lines and dots) of the map I>d 
(see Sec. IIV A 2[) and of the corresponding integrable system 
(thin red lines) of order K — 15. Here ffrcg accurately resem- 
bles the regular dynamics of U. (b) Orbits of the standard 
map (see Sec. IIVB|) for k — 2.9 (thick gray lines and dots) 
and of the corresponding integrable system (thin red lines) 
of order K = 7. Here -ffreg does not accurately resemble the 
regular dynamics of H. 



3. Quality of the prediction 

An important question is whether the direct tunnel- 
ing rates obtained using Eq. ([7|) depend on the actual 
choice of -ffrog and how these results converge in depen- 
dence of the order K of its perturbation series. There 
are two main problems: First the classical expansion of 
the Lie transformation is asymptotically divergent [63| . 
which means that from some order K on the series fails 
in reproducing the dynamics of the mixed system inside 
the regular region, see Fig. [2][b) . Second, for the quan- 
tization of ffreg its behavior in the vicinity of the last 
surviving KAM torus must be smoothly continued be- 
yond the regular island of U . Large fluctuations of i/rcg 
in this region, as appear for the method based on the fre- 
quency map analysis for large if, make the use of Eq. ([7]) 
impossible, see Fig.[3Ib). 




-0.5 0.0 q 0.5 -0.5 0.0 q 0.5 



FIG. 3. (Color online) Application of the method using the 
frequency map analysis. We show orbits of the standard map 
(see Sec. IIVB|) for k = 2.9 (thick gray lines and dots) and 
the corresponding integrable system (thin red lines) of order 
(a) if = 2 and (b) K = 10. While in (b) Jfreg resembles 
the regular island of U with higher accuracy than in (a), the 
extrapolation of -ffrcg beyond the island strongly oscillates. 
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Ideally one would like to use classical measures, which 
describe the deviations of the regular system i?rcg from 
the originally mixed one, to predict the error of Eq. ([7]) 
for the tunneling rates. However, these classical measures 
can only account for the deviations within the regular re- 
gion but not for the quality of the continuation of i?rcg 
beyond the regular island oiU . It remains an open ques- 
tion how to obtain a direct connection between the error 
on the classical side and the one for the tunneling rates. 

Nevertheless, the convergence of the integrable approx- 
imation can be studied by considering the tunneling rates 
determined with Eq. ([7|) under variation of the perturba- 
tion order K. For the example map I'd (introduced in 
Sec. IIV A 21) we find convergence up to the maximal con- 
sidered order. Fig. IH^a), and later use K = 10 for the 
comparison of Eq. ([7]) and numerical rates. For the stan- 
dard map at k = 2.9 the rates diverge rather quickly. 
Fig. Sfb), and we use K = 2. 

In general, different classical methods are applicable 
to determine a fictitious integrable system [/reg which 
leads to an accurate prediction of tunneling rates with 
Eq. (O. Hence, the determination of J/rog is not unique. 
The quality of an integrable system can be estimated a 
posteriori by comparison of the predicted tunneling rates 
with numerical rates. 



4- Application to billiards 

There are only few integrable two-dimensional billiard 
systems such as the circular, the rectangular, and the 
elliptical billiard. We use such integrable systems for 
various applications, as discussed in Sec. |V] for the mush- 
room and the annular billiard as well as for wires in a 
magnetic field and the annular microcavity. A general 
procedure to obtain i?rcg for arbitrary billiards is still 
under development. 
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C. Semiclassical evaluation 

In this section we semiclassically evaluate the direct 
rcgular-to-chaotic tunneling rates for systems where the 
fictitious integrable system is of the form H,-cg{q,p) — 
+ W{q). We consider one-dimensional kicked sys- 
tems 



H{q,p,t)=T{p) + V{q)Y,rS{t 



(18) 



which are the simplest Hamiltonian systems showing a 
mixed phase-space structure. They are described by the 
kinetic energy T{p) and the potential V{q) which is ap- 
plied once per kick period r = 1. The classical dynamics 
of a kicked system is given by its stroboscopic mapping, 
e.g., evaluated just after each kick 



Qn+l 
Pn+1 



qn 

Pn 



T'iPn 
' V'iqn 



(19) 



It maps the phase-space coordinates after the nth kick 
to those after the (n -|- l)th kick. The corresponding 
quantum map over one kicking period is given hy U = 
UvUt with 



Ut 



(20) 
(21) 



We consider a compact phase space with periodic bound- 
ary conditions for q € [-1/2, 1/2] and p S [-1/2, 1/2]. 

For an analytical evaluation of Eq. Q , which predicts 
the direct regular-to-chaotic tunneling rates, we approx- 
imate the fictitious integrable system U^eg by a kicked 
system, [/j-eg — Uyllr or C/rcg — UyUf with 



(22) 
(23) 



FIG. 4. (Color online) Tunneling rates 7 determined with 
Eq. ([7|, normalized by the numerical value 7num for m — 0, 
heS ~ 1/32 vs order K of //rog- In (a) we choose the system 
I'd (see Sec. IIV A2p and use the Lie transformation for the 

determination of Hrag, while in (b) the standard map (see ^rcg replaced by f/rog, giving 



Here the functions V{q) and T{p) are a low order Tay- 
lor expansion of V{q) and T{p), respectively, around the 
center of the regular island. Note, that the classical dy- 
namics corresponding to [/rcg is typically not completely 
regular. Still the following evaluation is applicable if U^cg 
has the properties: (i) Within the regular island it has an 
almost identical classical dynamics as U , including non- 
linear resonances and small embedded chaotic regions, 
(ii) It shows predominantly regular dynamics for a suf- 
ficiently wide region beyond the border of the regular 
island of U. 

We now give a semiclassical evaluation of Eq. ^ as- 
suming that both properties (i) and (n) are fulfilled. We 
consider the first case t/reg = Uyllr- As the dynamics of 
U and C/rcg are almost identical within the regular island 
of U, the approximate result, Eq. ([8]), can be applied with 



Sec. IIVB|I at K = 2.9 is considered for which we determine 
i^rcg using the method based on the frequency map analysis. 



\iu-UyUT)mg)f 

1)C/^(7t|C 



(24) 
(25) 
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We now use that \'ip'^^^ which is an eigenstate of the ex- 
act Uj-ae and i?reg, IS an approximate eigenstate of Uy-Ux, 



We obtain 



leading to C/^C/rlV-r^g) « e''^-«|V'r^g: 

im^\\{Uvui-t)w:i^)r. 

In position representation this reads 

( ^V{qk) 



N-l 



In 



2E |V'"g(ft)| 



fc=0 



1 



cff 



(26) 



(27) 



where AT/(g) :=F(g)-y(g) and ^k/N -1/2. In the 
semiclassical limit the sum in Eq. (j27l) can be replaced 
by an integral over the position space 



1/2 



-1/2 



cosf^) 



(28) 



Here, for the normalization J dq\il)'^J^q)\'^ = 1 holds, 
while previously |V'reg(9fc)P = 1 fulfilled. Note, 
that for the second case, where JJrcg ~ UyUf is used in 
Eq. (IS]), a similar result can be obtained in momentum 
representation. 



1/2 

2 J dp|V,;'g(p)| 

-1/2 



/AT(P)\ 

V fi-off / 



(29) 



with AT{p) := r(p) - T(p). 

We now use a WKB expression for the regular states 
[■f/;™ ). For simplicity we restrict to the case 



2 



W{q) 



leading to 



2^b(g)| 



exp 



V 



(30) 



(31) 



which is valid for q > q^. Here q^ is the right classi- 
cal turning point of the mth quantizing torus, w is the 
oscillation frequency, and p{q) 



I2{E^^^-W{q)). The 

eigenstates V'TOg('J') decay exponentially beyond the clas- 
sical turning point q^. The difference of the potential 
energies lS.V{q) approximately vanishes within the regu- 
lar region and increases beyond its border to the chaotic 
sea. Hence, the most important contribution in Eq. (|28|) 
arises near the left or the right border, q[ or qj, of the 
regular island. For q > ql we rewrite the regular states 

/ 

V'reg ilb) exp 

\ '"cn 



V'" (9) 



^yb(.')id.' 



^(32) 



V 

{ql)exp(-^{q-qlMql)\ 

V "cff 



(33) 



where in the last step we use p{q) ~ p{ql) in the vicinity 
of the border. 

In order to evaluate Eq. ([SS)) we split the integration 
interval into two parts, such that 7™ = 7^ + 7m j corre- 
sponding to the contributions from the left and the right. 
For simplicity we now approximate AV{q) by a piecewise 
linear function, 



AV{q) 



,qL<q< ql 

Cb{q-ql) , q> ql 



(34) 



with a constant Ch- With this we find 



7;„«2MV'.^g(<7DI' 



TT 

where x = (g - ql)/hcs, 



-2a:b(9nl 



[1 — cos(c6a;)] dx(35) 



(36) 



(1/2 - ql)/Ks. and 



[1 — cos(cf,a;)] dec. 



(37) 



In the semiclassical limit Xmax — ^ 00 and for fixed 
quantum number m the integral / becomes an hcft 
independent constant. The tunneling rate 7^ is propor- 
tional to the square of the modulus of the regular wave 
function at the right border q^ of the regular island. With 



Eq. (|3T|) we obtain 



luhcs 



2^^p{ql 



■ exp 



V 



2 

fi^cff 



b(9')|dg' 



(38) 



A similar equation holds for 7^. Note, that the same ex- 
ponent is obtained when considering the one-dimensional 
tunneling problem through an energy barrier in between 
the right turning point q^-^ and the right border of the 
island g^. 

As an example for the explicit evaluation of Eq. ([55)1 
we consider the harmonic oscillator Hreg{q,p) = p'^/2 + 
uj'^q'^/2, where oj denotes the oscillation frequency and 
gives the ratio of the two half axes of the elliptic invariant 
tori. Its classical turning points q^ = ±-y2^^/w, the 
eigenenergies = hca^{m -f 1/2), and the momentum 
\/2Era — q^oj'^ are explicitly given. Using these 
(|38)) and 7™ = 27^ we obtain 



expressions in Eq 



c — — exp 



9 A 

^^rcg 

hcs 



I3„ 



In 



(39) 

as the semiclassical prediction for the tunneling rate of 
the mth regular state, where ^^og is the area of the reg- 
ul ar islan d, a™ = (m + l/2){Arcg/hcs)~^ , and /3m = 
- am- The exponent in Eq. ([39]) was also derived in 
[2^ using complex-time path integrals. The prefac- 



Ref. 

tor 



(40) 
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can be estimated semiclassically by solving the integral, 
Eq. ([571) . for Xmax oo. For a fixed classical torus of 
energy E one obtains 



1 



With this prefactor the prediction Eq. gives excel- 
lent agreement with numerically determined rates over 
10 orders of magnitude in 7, see Fig. [TOl c). For a fixed 
quantum number m in the semiclassical limit the en- 
ergy Em approaches zero such that one can approximate 
b(9b)l ~ ^'i- which does not depend on /loff. 

Let us make the following remarks concerning Eq. p9p : 
The only information about this non-generic island with 
constant rotation number is Aiog/ft-cft as in Ref. .28j. In 
contrast to Eq. ([T]) it does not require further quantum 
information such as the quantum map U . While the term 
in square brackets semiclassically approaches one, it is 
relevant for large /loff. In contrast to Eq. (|28p . where the 
chaotic properties are contained in the difference AV{q), 
they now appear in the prefactor c via the linear approx- 
imation of this difference. 

In the semiclassical limit the tunneling rates predicted 
by Eq. p9|) decrease exponentially. For hctf — > one has 
a™ ^ and ^ 1, such that 7 ~ e"2'^«='5/'ioft, xhis 
reproduces the qualitative prediction obtained in Ref. @ . 
The non-universal constant in the exponent is 2 which 
is comparable to the prefactor 3 — ln4 1.61 derived in 
Refs. [28,, S]. We find that our result gives more accurate 
agreement to numerical rates, as will be shown in Sec. lIVI 
Still, a semiclassical evaluation of Eq. ([7]) for a fictitious 
integrable system of a more general form than Eq. ([30]) 
has to be developed. 



|65| with geometric mean [2l| 



(A(/3„i) = 



(42) 



(41) where describes the effective coupling of the mth 
regular state to the chaotic sea and r is the period of 
the driving. Note, that the factor y/Nch arises due to the 
different convention in Ref. [2l|, Eq. (1.27)], where the 
regular state is coupled to one chaotic state only. 

We now show that this average tunneling splitting 
(Aifm) in systems with symmetry- related regular regions 
is identical to the tunneling rate 7,^ from one regu- 
lar region to the chaotic sea: We start from Eq. (jA2p . 
7„i = Nch{\vch,m\'^) , and use the relation of the di- 
mensionless coupling matrix elements Vch m , defined in 

Eq. ®, to v:^ 



{\Vch- 



Together with Eq. ([42]) this leads to 



(43) 



7m = {Aipm)- 



(44) 

This result was previously employed in Fig. 3 in Ref. [s^ , 
where numerical splittings Aipm are used and the predic- 
tion is for tunneling rates 7™. 

Fig. [5] illustrates the strong fluctuations of the split- 
tings A</?„i (squares) in contrast to the smooth behavior 
of the tunneling rates 7™ (dots, lines). As predicted by 
Eq. (|44| one can see in the figure that the splittings fluc- 
tuate around the tunneling rates as a function of 1/ h^s- 

This demonstrates that for a quantitative verification 
of a theory on regular-to-chaotic tunneling the tunneling 
rates allow for a more precise comparison than tunneling 
splittings. 



D. Relation to chaos-assisted tunneling 

In Ref. j9j Tomsovic and UUmo studied dynamical tun- 
neling in systems with two symmetry-related regular is- 
lands surrounded by a chaotic region in phase space. 
They considered the quasi-energy splittings Aipm be- 
tween the symmetric and antisymmetric regular states 
on the mth quantizing tori of both islands. These tun- 
neling splittings are drastically enhanced by the presence 
of chaos, i.e. chaotic states assist the tunneling process 
compared to the case of a system with integrable dynam- 
ics between the regular islands. The two-step process, 
which couples the regular torus from one island to the 
chaotic sea and from the chaotic sea to the symmetry- 
related torus of the other island, dominates the direct 
coupling of the two regular tori. 

The tunneling splittings A(^m show fluctuations over 
several orders of magnitude under variation of external 
parameters [1, H, [s^l • These fluctuations originate from 
the varying distance of the regular doublet to the chaotic 
states and their varying coupling. According to a random 
matrix model, the splittings follow a Cauchy distribution 




10 30 50 l/h^s 70 



FIG. 5. (Color online) Dynamical tunneling rates 7 (dots), 
quasi-energy splittings Aip (squares), and the prediction of 
Eq. ((Tj) (lines) vs l//icft for the regular states m < 2. We 
use the system I'd (see Sec. IIV A2|) . The insets show the 
phase space with two symmetry-related regular regions (used 
to determine Ai/j) and the phase space with one regular region 
(used to determine 7). 
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III. NUMERICAL DETERMINATION OF 
TUNNELING RATES 

To test the theoretical prediction derived in Sec. Illlwe 
compare its results to numerical rates in Sees. IIVI and 
fVl In this section we present three alternative methods 
to numerically compute tunneling rates: (A) opening the 
system, (B) time evolution of regular states, and (C) eval- 
uating avoided crossings. Figure [5] shows a comparison of 
the tunneling rates obtained by these three methods for a 
quantum map. We find excellent agreement between the 
first two methods while the last approach shows small 
deviations. 



A. Opening the system 

The structure of the considered phase space, with one 
regular island surrounded by the chaotic sea, allows for 
the determination of tunneling rates by introducing ab- 
sorption somewhere in the chaotic region of phase space. 
For quantum maps this can be realized, e.g., by using a 
non- unitary open quantum map (66l [67j 



U° = PUP, 



(45) 



where P is a projection operator onto the complement of 
the absorbing region. An example is given by a sum of 
projectors on position eigenstates. 



p = J2\i)(i\^ 



(46) 



where the regular island is located well inside the interval 

[qi,Qr]- 

While the eigenvalues of U are located on the unit 
circle the eigenvalues of U° are inside the unit circle as 



10" 
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H B^"h 

B " i ■ 



U° is sub-unitary, see Fig. [T] The eigenequation of U° 
reads 



with eigenvalues 



(47) 



(48) 



The decay rate is characterized by the imaginary part of 
the quasi-energics in Eq. and one has 



7m = -21og|z„ 



2(1 



J). 



(49) 



In order to obtain the open map U° practically, we 
quantize the classical map on the cylinder {q,p) € 
[— oo, oo] X [— 1/2, 1/2] with the periodically extended po- 
tential V{q). This leads to an infinite dimensional uni- 
tary matrix U in position representation (68j and we find 
U° with Eq. using the projector P given by Eq. (|46p . 
After the diagonalization of U° we identify the eigenval- 
ues Zm of U° close to the unit circle, which correspond 
to the quasi-bound regular states, and use Eq. (|49|) to 
determine the tunneling rates, see Fig. [S] (dots) for an 
example. 

If the chaotic region does not contain partial barriers 
and shows no dynamical localization, it is justified to as- 
sume that the probability of escaping the regular island is 
equal to the probability of leaving through the absorbing 
regions located in the chaotic sea. Then, the location of 
the absorbing regions in the chaotic part of phase space 
has no effect on the decay rates. 

In generic systems, however, partial barriers will ap- 
pear in the chaotic region of phase space. The addi- 
tional transition through these structures further limits 
the quantum transport such that the calculated decay 
through the absorbing region occurs slower than the de- 
cay from the regular island to the neighboring chaotic 
sea. Similarly, dynamical localization in the chaotic re- 
gion may slow down the decay. The quantitative infiu- 
ence of partial barriers and dynamical localization on the 
regular-to-chaotic tunneling rates is an open problem for 
future studies. If necessary we will suppress their infiu- 
ence by moving the absorbing regions closer to the regular 
island. 



(c) 




1 10" 


-H- ++'+++ + 


l-\z\ 


+ 

* 


10"^ 


* 

* 


10"'' 





10 



20 



30 



40 



50 60 



Re 2 1 



Re 2 







FIG. 6. (Color online) Comparison of tunneling rates 7 ob- 
tained by opening the system (dots), time evolution (crosses), 
and the evaluation of avoided crossings (squares) vs 1/ftcff- 
We use the map ©ho (see Sec. IIV A If) and consider the reg- 
ular states m < 3. The inset shows the phase space of the 
system. 



FIG. 7. (Color online) (a) The eigenvalues of the unitary 
quantum map ©ho (see Sec. IIV A 1|) with h^s ~ 1/14 are lo- 
cated on the unit circle. Regular states are marked by red 
stars, chaotic states by blue plus symbols, (b) The eigenval- 
ues of the open system U° are located inside the unit circle, 
(c) The distance 1 — 1^1 is shown on a logarithmic scale. 
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B. Time evolution 

A simple method to obtain a numerical prediction of 
the tunneling rates for quantum maps is given by the 
time-evolution of a purely regular state IV'reg) with a non- 
unitary operator U" = PUP. Here P projects onto a 
region in phase space which includes the regular island. 
We consider 

for t G N, which describes the probability of the time- 
evolved regular state in the regular island at time t. At 
each time step some probability of IV'reg) is absorbed 
in the chaotic region due to the openness of the quan- 
tum map U°. Consequently, Wm{t) decays exponentially, 
Wm{t) ~ e"'''"*, and the tunneling rates 7™ can be de- 
termined by a fit of the numerical data. If IV'reg) con- 
tains admixtures from lower excited regular states (with 
smaller tunneling rates) their decay dominates at times 
t 3> l/7m- If it contains admixtures from higher excited 
regular states (with larger tunneling rate) their decay will 
be seen at small times, see Fig. [Sj The computed tun- 
neling rates are in excellent agreement with the results 
obtained by opening the system, see Fig. [S] (crosses). This 
method works best for regular states IV-'reg) which resem- 
ble the corresponding eigenstates of the mixed system U 
with high accuracy. It is particularly useful for Hilbert 
spaces of large dimension N where diagonalizing the ma- 
trix U° would numerically be very time-consuming. 



C. Evaluation of avoided crossings 

The third method calculates the tunneling rate of a 
regular state directly from the spectrum of the system. 
For quantum maps we determine the quasi-energies ip un- 
der variation of a parameter of the system which leaves 




t 100000 t 100000 

FIG. 8. (Color online) Time evolution of a purely regular state 
for the system ©ho (see Sec. IIVAII) with /i^ff = 1/14. We 
present W{t) vs t obtained by Eq. (|50|) on a semi-logarithmic 
scale, (a) For the regular ground state the tunneling rate to 
the chaotic sea 70 is determined by the slope of the numerical 
data, (b) For the first excited state m — 1 the slope for small 
times t determines the tunneling rate 71, while for larger times 
the decay of the ground state 70 with a smaller slope is found. 



the classical dynamics invariant, such as the Bloch phase 
6g or 9p. These phases specify the periodicity conditions 
on the torus and can be incorporated in the quantiza- 
tion of the map f6^. Under variation of such a pa- 
rameter the quasi-energy of the considered regular state 
(pm shows avoided crossings with quasi-energies tpch of 
chaotic states, see Fig. [S] These avoided crossings have 
widths Aipch,m- According to degenerate perturbation 
theory they are related to the matrix elements Wch,m by 
'-^Vch,m = 2wch,m and fluctuate depending on the in- 
volved chaotic state. The tunneling rate follows from the 
dimensionless version of Fermi's golden rule, Eq. (IA2p . 

7m = ^^{\Aipch,m\'^)ch, (51) 

where A'ch is the number of chaotic states. Note, that 
the two methods discussed in Sees . IIII Al and IIITbI deter- 
mine the tunneling rates 7 for fixed Bloch phases. Under 
variation of 9q or 9p we observe numerically for these 
methods that the rates vary by at most a factor of 2, 
while Eq. ([FT]) gives an average rate. 

The quality of this prediction depends on the number 
of avoided crossings entering the average in Eq. ([51]) . If 
only a few avoided crossings are included, the statisti- 
cal error of the result is large. Note, that for quantum 
maps with a mean drift in the chaotic sea, such as the sys- 
tems V (introduced in Sec. lIV A"T|) . under variation of the 
Bloch phases 6 several chaotic states will show avoided 
crossings with each regular state fE§\ ■ The results of this 
method are presented in Fig. |6] (squares) for an example 
system. We typically obtain tunneling rates which are 
smaller than the results of the other two methods. While 
we have no explanation of this behavior, the deviation 
is smaller than a factor of two, which is sufficient for a 
comparison to theoretical predictions. 

Also for billiards tunneling rates can be computed by 
this method. We determine the spectrum under varia- 
tion of parts of the billiard boundary which leaves the 
classically regular dynamics unchanged but affects the 




FIG. 9. (Color online) Quasi-energies of the quantum map U 
for the system ©ho (see Sec. IIV A 1[) vs Bloch phase 9q with 
ftcff = 1/14 and 6p = 0. Regular states (red solid lines) show 
avoided crossings with chaotic eigenstates (blue dashed lines). 
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chaotic dynamics. Quantum mechanically, the eigenener- 
gies of the regular states remain almost unaffected while 
the eigenenergies of the chaotic states vary strongly, due 
to the changing density of chaotic states. Hence, avoided 
crossings of widths AEch,mn between the considered reg- 
ular and the chaotic states appear. The tunneling rate is 
given by Fermi's golden rule, Eq. (ITil) . 



,mn I Pch) 

Y 4 



16M 



(|Ai?ch,mnP-4ch) 



(52) 

where we average the product of all numerically deter- 
mined widths Ai?ch,iTm = 214h, mn ^ud the corresponding 
density of chaotic states pch, which we approximate by 
its leading Weyl term, see Sec. Ill A 31 Note, that in gen- 
eral it can be difficult to deform a part of the billiard 
boundary such that the regular dynamics is unchanged 
while still the numerical methods for the determination 
of eigenvalues in billiard systems are applicable. 



IV. APPLICATION TO QUANTUM MAPS 

In the following we will apply the fictitious integrable 
system approach, derived in Sec. [TTl to the prediction 
of direct regular-to-chaotic tunneling rates in the case 
of quantum maps and compare the results to numerical 
rates for different example systems. 



A. Designed maps V 

Our aim is to introduce kicked systems which can be 
designed such that their phase space shows one regular 
island embedded in the chaotic sea, with very small non- 
linear resonance chains within the regular island, a neg- 
ligible hierarchical region, and without relevant partial 
barriers in the chaotic component. For such a system 
it is possible to study the direct regular-to-chaotic tun- 
neling process without additional effects caused by these 
structures. 

To this end we define the family of maps V, according 
to Eq. , with an appropriate choice of the functions 
T'{p) and V'{q) [H [al HI, in, m. For this we first 
introduce 



t'ip) 
v'{q) 



i - (1 - 2p) for 
| + (l-2p) for 

-rq + Rq^ for 



~i < P < 



(53) 
(54) 



with < r < 2 and i? > 0. This gives a regular island 
around [q^p) ~ (0, 1/4). Considering periodic boundary 
conditions the functions t' {p) and v'{q) show discontinu- 
ities at p = 0, ±1/2 and q = ±1/2, respectively. In order 
to avoid these discontinuities we smooth the periodically 
extended functions v'{q) and t'{q) with a Gaussian, 



G{z) 



1 



V2 



: exp 



z 

2^ 



resulting in analytic functions 

T'{p)^ j Azt'{z)G{p~ z), (56) 
V'{q)^ [ dzv'{z)G{q~z), (57) 



(55) 



which are periodic with respect to the phase-space unit 
cell. With this we obtain the maps T) depending on the 
parameters r, R, and the smoothing strength e. The 
smoothing e determines the size of the hierarchical region 
at the border of the regular island. Tuning the parame- 
ters r and R one can find situations, where all nonlinear 
resonance chains inside the regular island are small. 



1. Map T)ho with harmonic oscillator-like island 

For R — Q both functions v'{q) and t'{p) are linear in 
q and p, respectively. In this case we find a harmonic 
oscillator-like regular island with elliptic invariant tori 
and constant rotation number. We choose the parame- 
ters r = 0.46, i? = 0, e = 0.005 and label the resulting 
map by I?ho- Its phase-space is shown in the insets of 
Fig. 1101 Numerically, we determine tunneling rates by 
introducing absorbing regions at |g| > 1/2, as described 
in Sec. IIII Al In order to apply the fictitious integrable 
system approach we use the Hamiltonian of a harmonic 
oscillator as iJrcg- It is squeezed and tilted according 
to the linearized dynamics in the vicinity of the stable 
fixed point located at the center of the regular island. 
Its eigenfunctions IV'reg) analytically known, see Ap- 
pendix |B] 

Figure llOf a) shows the numerically evaluated predic- 
tion of Eq. ([7]) compared to numerical tunneling rates. 
We find excellent agreement over more than 10 orders of 
magnitude in 7. In the regime of large tunneling rates 
small deviations occur which can be attributed to the in- 
fluence of the chaotic sea on the regular states: These 
states are located on quantizing tori close to the border 
of the regular island and are affected by the regular-to- 
chaotic transition region. However, the deviations in this 
regime are smaller than a factor of two. 

Figure fTOl b) shows the results of Eq. ([ST)) , which are 
obtained by approximating Urcg by a kicked system, 

f^reg ~ UyUr, again using the analytically given IV'rcg)- 
These results are still in excellent agreement with the nu- 
merical rates (solid lines). In Eq. (pS)) the sum over the 
positions is replaced by an integral, which explains the 
small deviations to the results of Eq. (l27l) , see Fig. ITOl' b) 
(dashed lines). These deviations vanish in the semiclas- 
sical limit. 

Finally, in Fig. fTUTc) we compare the results of the 
semiclassical prediction, Eq. to the numerical rates. 
Due to the approximations performed in the derivation of 
this formula stronger deviations are visible in the regime 
of large tunneling rates while the agreement in the semi- 
classical regime is still excellent. 
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FIG. 10. (Color online) Numerical tunneling rates (dots) for 
m < 8 for the map I>ho with a harmonic oscillator-like island, 
(a) Comparison with Eq. (O (lines). The insets show Husimi 
functions of the regular states for the quantum numbers m = 
and m = 5 at 1/hcS = 30 and the classical phase space 
of the system, (b) Comparison with Eqs. H27p (solid lines) 
and (|28p (dashed lines), (c) Comparison with Eq. (I39|) (solid 
lines). The prediction of Refs. 28, 57], Eq. ((58)l . for m = 
with a fitted prefactor is shown (dotted line). 



In Refs. [28|, |57| a prediction was derived for the tun- 
neling rate of the regular ground state. 



70 



r(Q, 4q) 
' r(a,o) 



where V is the incomplete gamma function, a = 
^rcg/^off, and c is a constant. Equation (|58|) can be ap- 
proximated semiclassically [2l|, a — > oo, leading to 



1 



70 cx 



-Q(3-ln4) 



(59) 



Figure [TUTc) shows the comparison of Ec^. (1551) (dotted 
line) to the numerical rates for the map I?ho- Especially 
in the semiclassical regime strong deviations are visible. 
The factor 2 which appears in the exponent of Eq. ((39)) 
is more accurate than the factor 3 — In 4 in Eq. (IST 



2. Map Vd with deformed island 

In generic systems the regular island has a non-elliptic 
shape and the rotation number of regular tori changes 
from the center of the regular region to its border with 
the chaotic sea. Such a situation can be achieved for the 
family of maps V with the parameter i? ^ 0. For most 
combinations of the parameters r and R resonance struc- 
tures appear inside the regular island. They limit the 
/loff-regime in which the direct regular-to-chaotic tunnel- 
ing process dominates. Hence, we choose a situation in 
which the nonlinear resonances are small such that their 
influence on the tunneling process is expected only at 
large l/heg. For this we use r = 0.26, R = 0.4, e = 0.005 
and label the resulting map with a deformed island by 
I?d, see the inset in Fig. [TT] for its phase space. 

We determine the fictitious integrable system H^eg by 
means of the Lie-transformation method described in 
Sec. Ill Bl It is then quantized and its eigenfunctions are 
determined numerically. Figure [TT] shows a comparison 
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FIG. 11. (Color online) Dynamical tunneling rates from a 
regular island to the chaotic sea for the map Vd- Numerical 
rates (dots) and prediction following from Eq. (O (lines) vs 
1/heH for quantum numbers m < 5. The insets show Husimi 
representations of the regular states m = and m = 5 at 
1/hes = 50. The prediction of Refs. [H, [53|, Eq. dSH), for 
m — with a fitted prefactor is shown (dotted line). 
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of the numerically evaluated prediction of Eq. ([7]) (solid 
lines) to numerical tunneling rates (dots) yielding excel- 
lent agreement for 7 > 10^^^. For smaller values of 7 de- 
viations occur due to resonance-assisted tunneling which 
is caused by a small 10:1 resonance chain. Similar to the 
case of the harmonic oscillator-like island the fictitious 
integrable system Ureg can be approximated by a kicked 
system Urcg « UyllT using V{q) = -rq^/2 + 
Hence, Eqs. (fST)) and (|28l) can be evaluated giving sim- 
ilarly good agreement (not shown). The prediction of 
Eq. ((55)) [28, 57] (dotted fine) shows large deviations to 
the numerical rates. 



3. Map X>™c with weakly chaotic dynamics 

In Ref. 30] the dynamical tunneling process from one 
regular island to the chaotic sea in a system with weakly 
chaotic dynamics was investigated. Here, tunneling can 
occur to regions in phase space which are far away from 
the border of the regular island with the chaotic sea. 
We show that also in this situation, which we believe to 
be non-generic, the fictitious integrable system approach 
can be applied. Its results will be compared to the WKB 
prediction of Ref. Q • 

A system with weakly chaotic dynamics can be mod- 
eled by the example systems V. We choose r = 0.05, 
R = 0.1, and e = 0.005, consider the extended phase 
space {q,p) S [—1, 1] x [—1/2, 1/2], and label the result- 
ing map by 2?wc- Its phase space is shown in the upper 
inset of Fig. [ll 

In order to apply the fictitious integrable system ap- 
proach we use the Lie-transformation method, as de- 
scribed in Sec. Ill Bl to obtain the fictitious integrable 
system iJiog- As the system P^c is only weakly driven, 
due to the small parameters r and R, it is sufficient to 
consider the zeroth order of the Lie expansion which has 
no mixed terms containing q and p simultaneously. The 
resulting integrable approximation 

Hrog{q,p)^^ + W{q), (60) 

describes the dynamics in a potential W{q) = u'^q^l2 - 
Rq'^/6 with oj = ^ r/2, see the lower inset in Fig. [12] 

In Ref. [131 it was shown that for such a weakly chaotic 
system regular-to-chaotic tunneling rates can be pre- 
dicted by one-dimensional tunneling under the energy 
barrier of the potential W{q). For the tunneling rates 
one finds 

7™«^exp[-|- Jb(g,i?™)|dJ (61) 



where p(q, E"^^) = ^'^^"4 ^ ^^^^^ + ^9^3, q'm denotes 
the right classical turning point inside the potential well, 
gj^ is the turning point outside the potential well, and 
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FIG. 12. (Color online) Dynamical tunneling rates from a 
regular island to the chaotic sea for the weakly chaotic system 
fwc- We compare numerical rates (dots) and the prediction 
following from Eq. ^ (solid lines) and Eq. (|6l]) [301 (dashed 
lines) vs l/h^s for m = and m = 1. The upper inset shows 
the phase space of the system, where the absorbing regions 
are indicated in gray. In the lower inset the approximate one- 
dimensional potential W{q) used in Eq. (I6ip is presented with 
the inner (outer) turning point ql^ (qm) at energy i?™g. 



LOm is the oscillation period on the mth quantizing torus. 
The eigenenergies can be calculated using the Bohr- 
Sommerfeld quantization § p{q, E^^) dq = hef{{m+ 1/2). 
For the system 2?wc the right turning point 5^ is located 
far away from the regular island. Tunneling occurs to 
the region with q > deep inside the weakly chaotic 
sea and not to the neighborhood of the regular island, as 
for the other examples considered in this paper. 

In Fig. [T^j we compare the numerically evaluated pre- 
diction of Eq. ([7|) (solid lines) to the result of Eq. ([6T|) 
(dashed lines) and numerical rates (dots), which are de- 
termined by absorbing regions at \q\ > 1. We find good 
agreement. 

Note, that for the weakly chaotic system P^c the 
purely regular states IV'reg) ^reg show the correct tun- 
neling tails far beyond the regular island including the 
outer turning point q'^. Moreover, the semiclassical eval- 
uation of Eq. ([8]), presented in Sec. Ill C[ can be per- 
formed. This leads to Eq. ([55)) which has the same expo- 
nential term as Eq. ()6ip but a different prefactor. 

We want to emphasize that generically regular-to- 
chaotic tunnehng cannot be described by Eq. ([CT)) . as the 
integrable approximation Hj-cg is not of the form Eq. ()60p . 

B. Standard map 

The paradigmatic model of an area preserving map 
is the standard map [zO], defined by Eq. ([T9)) with the 
functions 

T'{p)=p, (62) 
y'(g) = ^sin(2^g). (63) 
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For K between 2.5 and 3.0 one has a large generic regu- 
lar island with a relatively small hierarchical region sur- 
rounded by a 4:1 resonance chain, see the inset in Fig.lTSl 

When determining tunneling rates numerically by in- 
troducing absorbing regions at \q\ > 1/2 we find strong 
fluctuations as a function of ft-cff, presumably caused by 
partial barriers. Using absorption at \q\ > 1/4, which is 
closer to the island, we find smoothly decaying tunneling 
rates (dots in Fig. [13]). 

Evaluating Eq. ([7|) for k = 2.9 gives reasonable agree- 
ment with these numerical rates with deviations up to a 
factor of 5, see Fig. [T3] (solid lines). Here we determine 
TJreg using the method based on the frequency map anal- 
ysis as the Lie transformation is not able to reproduce 
the dynamics within the regular island of U, see Sec. Ill Bl 
With increasing order K of the expansion series of //reg 
the tunneling rates following from Eq. ([7]) diverge, see 
Fig.m^b). Hence, for the predictions in Fig. [13] we choose 
K = 2 which is the largest order before the divergence 
starts to set in. Note, that at such small order K the 
accuracy of TJ^eg within the regular region of U is infe- 
rior compared to the examples discussed before. Hence, 
in Eq. ([7]) the state U\tp^^) has small contributions of 
other purely regular states IV'reg) regular island. 

These contributions are compensated by the application 
of the projector Pch- However, this projector depends on 
the number of regular states A^rog, which grows in the 
semiclassical limit. If A^rcg increases by one, Pj-cg sud- 
denly projects onto a larger region in phase space. This 
explains the steps of the theoretical prediction, Eq. ([7]), 
visible in Fig. 1131 How to improve -ffrcg and the projector 
Pch is an open question. 

C. Map Vra with a regular stripe 

Another designed kicked system was introduced in 
Refs. [13, [H, [Zll, Izll- Here the regular region consists 
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FIG. 13. (Color online) Tunneling rates for the standard map 
(k = 2.9) for m < 2 vs l/heS- Prediction of Eq. (lines) 
and numerical rates (dots), obtained using absorbing regions 
£^t |g| > 1/4 (gray-shaded area of the inset). 



of a stripe in phase space, see the inset in Fig. [Ml In our 
notation the mapping Pj-s, Eq. ()19|) . is specified by the 
functions 

V'{q) = -^{8TTaq + di-d2 

ZTT 

-|--[87rag — uj + di] tanh[b{8TTq — qd)] 
+ ^[-8TTaq + uj + d2] tanh[6(87rg + qa)]), (64) 
T'(p) = -^sin(2^p) (65) 

with parameters a — 5, b — 100, di = —24, d2 — —26, 
Lo = 1, qd — 5, and K — S. The kinetic energy T{p) is 
periodic with respect to the phase space unit cell. 

The resulting map 2?rs is similar to the system I?ho as 
it also destroys the integrable region by smoothly chang- 
ing the function V'{q) at |g| = qd/{8n). For \q\ < qd/ {8t:) 
the potential term is almost linear while it tends to the 
standard map for |g| > qd /{8Tr). The parameter b deter- 
mines the width of the transition region. In Ref. [7l| this 
map is used to study the evolution of a wave packet ini- 
tially started in the regular region by means of complex 
paths. We now predict direct regular-to-chaotic tunnel- 
ing rates with Eq. ([7])- The fictitious integrable system 
t/rcg is determined by continuing the dynamics within 
l^l < qd/{8n) to the whole phase space. It is given as a 
kicked system, Eq. ([T9]) . defined by the functions 

T'{p)^-^sin{2iTp). (67) 

When determining tunneling rates numerically using 
absorbing regions at \q\ > 1/2 we find strong fluctua- 
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FIG. 14. (Color online) Dynamical tunneling rates from a 
regular stripe to the chaotic sea for the map Drs- We com- 
pare numerical rates (dots) and the prediction following from 
Eq. ([7]) (lines) vs l//icff for the quantum numbers |m| < 4. 
The inset shows the phase space of the system. The numer- 
ical rates are obtained using absorbing regions at |g| > 1/4 
(gray-shaded area of the inset). 
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tions as a function of hes, similar to the standard map. 
Choosing \q\ > 1/4 for the opening, which is closer to 
the regular stripe, we find smoothly decaying tunneling 
rates (dots in Fig. I14p . Their comparison with the nu- 
merical evaluated prediction of Eq. ([7]) shows excellent 
agreement, see Fig. [14] (lines). 

Note, that due to the symmetry of the map there are al- 
ways two regular states with comparable tunneling rates 
except for the ground state m — 0. These two states 
are located symmetrically around the center of the regu- 
lar stripe. While the prediction, Eq. ([7]), is identical for 
both of these states, the numerical results differ slightly 
due to the different chaotic dynamics in the vicinity of 
the left and right borders of the regular region. 



V. APPLICATION TO BILLIARDS 

Billiards play a central role in both experimental and 
theoretical studies in quantum chaos. They are dynam- 
ical systems given by a point particle of mass M which 
moves with constant velocity inside a domain f2 e 
which we assume to be compact. The particle is elas- 
tically reflected at the boundary dfl such that the an- 
gle of incidence equals the angle of reflection. While 
there are only a few integrable and completely chaotic 
billiards, the majority shows a mixed phase space con- 
sisting of regions of regular and chaotic dynamics. Quan- 
tum mechanically, billiards are described by the time- 
independent Schrodinger equation (in units h ~ 2M ~ 1 
used in this section) 

- AV-r, {q) = Er,i>n{q), q^fl (68) 

with the Dirichlet boundary condition ipniq) = 0, q G 
dfl. In Eq. ([BS]) A denotes the Laplace operator in 
two dimensions. Equation is identical to the eigen- 
value problem of the two-dimensional Helmholtz equa- 
tion which for example describes electro-magnetic modes 
in a microwave cavity. This equivalence allows for the 
simulation of quantum billiards by experiments using mi- 
crowave cavities [73l - [77j . 

The state of a particle is described by a wave function 
tp{q) € L'^{^l) in position representation, where L^{il) 
is the Hilbert space of square integrable functions on ft. 
Due to the compactness of fl the eigenvalues {£"„} are 
discrete and can be ordered as < Ei < E2 < < ■ ■ ■ . 
The eigenfunctions can be chosen real and form an or- 
thonormal basis on L'^{il). In contrast to the case of 
quantum maps discussed in Sec. IIVI one gets infinitely 
many eigenvalues and eigenfunctions. There are only a 
few billiard systems, whose eigenfunctions are analyti- 
cally known, e.g., the rectangular, circular, and elliptical 
billiard. Usually an analytical solution of Eq. (|68)) is not 
possible. 

The determination of tunneling rates for two- 
dimensional billiard systems is of current interest. It is 
relevant, e.g., in the context of light emission in optical 



microcavities [47l - l49l . [TSf , flooding of regular states [i^ - 
lil ] and conductance properties of electrons in disordered 
wires with a magnetic field [H, [sgI . Previous theoreti- 
cal pred ictions of tunneling rates [31] or energy-splittings 
[H [i3 in billiards required additional free parameters. 

We apply the fictitious integrable system approach in 
order to determine direct regular-to-chaotic tunneling 
rates for billiards. For this we employ Eqs. (1131) and 
([T^ , where in the following we omit the tilde of the non- 
orthogonal chaotic states ipch and label the correspond- 
ing dimensionless matrix elements and tunneling rates 
by Well, mn and jmn, respectively. For the chaotic states 
entering in Eq. (jl3p we employ random wave models [3| 
such that the average in Eq. ([Ti)) becomes an ensemble 
average over the different realizations of the random wave 
model. From this we obtain explicit analytical predic- 
tions for the mushroom billiard [33| , the annular billiard, 
two-dimensional nanowires with one-sided surface disor- 
der, and optical microcavities ^34i] . 



A. Mushroom billiard 

We consider the desymmetrized mushroom billiard 
(79{ , see Fig. [TSla). characterized by the radius R of 
the quarter circular cap, the stem width a, and the 
stem height /. This billiard is of great current inter- 
est [31I, [51I, [80l - l82| due to its sharply separated regular 
and chaotic regions in phase space. The regular trajecto- 
ries show whispering-gallery motion and do not cross the 
small quarter circle of radius a. Each trajectory which 
crosses this curve is chaotic, see Fig. [WcV There is 
no hierarchical regular-to-chaotic transition region and 
there appear no resonance chains inside the regular is- 
land. Hence, for this billiard resonance-assisted tunneling 
does not occur and the direct regular-to-chaotic tunnel- 
ing process is relevant for all energies E. The application 
of the fictitious integrable system approach to the mush- 
room billiard leads to the explicit analytical formula (|83p , 
which was obtained in Ref. [33| , where it was successfully 
compared to experimental data. 



1. Derivation of tunneling rates 

In order to predict tunneling rates for the mush- 
room billiard we review the derivation [33] starting from 
Eqs. and First we construct a fictitious in- 

tegrable system iJreg, determine its eigenstates V'reg 
and find a model for the chaotic states ipch{q)- In the 
following analysis we set R = 1. A natural choice for 
the regular system 7?reg is the quarter-circle billiard. Its 
eigenfunctions are analytically known 

^/'"g (r, ip) = NranJ-miimnr) sin(TO.^), (69) 

in polar coordinates (r, Lp). They are characterized by the 
radial [n = 1,2,...) and the azimuthal (m = 2, 4, ... ) 
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where A depends on W via the Schrodinger equation as 

-\^ + W = EZn. (73) 

Since the regular eigenfunctions ip^^w ^^"^ their deriva- 
tives have to be continuous at y = we obtain 



A 



V'™gV(^,2/ = 0) 



(74) 



Evaluating Eq. (1131) for the coupling matrix elements one 
finds 





Vch.mn= lim Idx dyipch{x,y){-W)'ip'^^]y{x,y) (75) 

W-^oo J ' 








= — lim 



W 



dx / dy^eh(a;,y)— e^^5,^™gXx,0)(76) 



FIG. 15. (Color online) (a) Schematic picture of the mush- 
room billiard, with cap radius R, stem width a, and stem 
height /, showing the two coordinate systems used in the the- 
oretical derivation of the direct regular-to-chaotic tunneling 
rates, (b) Auxiliary billiard H^^. (c) Poincare section at the 
quarter-circle boundary (relative tangential momentum p vs 
arclength s) showing regular and chaotic regions with illus- 
trations of trajectories. 



In the last equation the term M^e'^^/A appears, which in 
the limit W ^ oo gives for y < 



W 



W 



,eJ^^^^v ^ 2<5(2/), (77) 



where we use Eq. ^ and that £:^„ remains bounded. 
For the coupling matrix elements, Eq. (|76p . we obtain 



m ; ^ ran 



quantum numbers. As we are considering the quarter- 
circle billiard m is allowed to take even values only. 
Here denotes the mth Bessel function, jVrm the nth 
root of J, 



/8/7r/Jm_i(jm„) accounts for the 
normalization, and Emn = Jmn ^"^^ the eigenenergies. 
Among the regular states of the quarter-circle billiard 
we will consider only those which concentrate on regular 
tori of the mushroom billiard with angular momentum 

Vmn — '^^l jran ^• 

We use the Hamiltonian H of the mushroom and iJrcg 
of the quarter-circle billiard in Eq. (fT3|) to determine the 
coupling between the regular and the chaotic states. An 
infinite potential difference H — Hj-cg ~ ~oo occurs, while 
V'rag = in the stem of the mushroom for y < 0. In 
order to avoid the undefined product {H ~ H^-cgj^p^g we 
introduce a finite potential at y < 0, see Fig. [TSTb). 



V{q) 



p2 + V{q) (70) 

for + y^ < 1, X, y > 
W for y < 0, < a; < 1 (71) 
oo otherwise 



and consider the limit W ~^ od in which the quarter- 
circle billiard is recovered. For finite W the regular eigen- 
functions V'rog wC*?) '^^ -^rog decay into the region y < 0. 
To describe this decay we make the following ansatz 



V'™gV(a:,y)=V^;™,v^(a;,y = 0)e 



(72) 



- / dxAhix,y^ 0)9,^""(a;,y = 0). (78) 



Due to the limiting process only an integration along the 
line y = remains which connects the quarter circle 
billiard to the stem of the mushroom. Equation ((78)) 
contains the derivative of the regular wave function per- 
pendicular to this line. The largest contribution of the 
integral is close to the corner of the mushroom at a; = a, 
as the derivative of the regular eigenfunctions dyip^g de- 
cays toward a; = 0. Inserting the regular states, Eq. ((M|) . 
one finds 



dxipch{x,y ^ 0) — Jm{jninX). (79) 

X 



In order to evaluate Eq. (j79p we use a random wave de- 
scription to model the chaotic states ipchiQ)- H has to re- 
spect the Dirichlet boundary conditions in the vicinity of 
the corner at a; = a. For this random wave model we use 
polar coordinates q = (p, i)) as introduced in Fig. [TSla) 
such that the corner of angle 37r/2 is located at (0,0). 
The Dirichlet boundary conditions at this corner are ac- 
counted for using (83| 



3 j:c.J.(Msin (80) 

S—1 



2s 



in which x — a = pcos(??) and y = psui('d). The co- 
efficients Cs are independent Gaussian random variables 
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with mean zero, (cg) = 0, and unit variance, {CgCt) = 5st- 
Equation ([50)1 fulfills the Schrodinger equation at en- 
ergy E — and the prefactor is chosen such that 
(|V'ch(/3, 'J)!^) = l/-4ch holds far away from the corner. 
Note, that we do not require these chaotic states to de- 
cay into the regular island, as Eq. ([79l) is an integral 
along a line of the billiard where the phase space is fully 
chaotic. Near the boundary, but far away from the cor- 
ner, (l^chP) recovers the behavior 1 — Jo{2k\x\) [s^-lssj. 
Inserting Eq. dHOl) into Eq. ([73) at energy E^n = Jmn 
and 1) = n, we obtain the coupling matrix elements 



10" 



Vch,i 



2s 

, Cs sm I TT 

3^ch ^ V 3 

S— 1 



10" 



10" 



10-1^ 



10" 




20 



40 



60 k 80 



J^{jmn[a- x]) — Jrn{jmnX)dx, (81) 
3 X 



where terms with s a multiple of 3 vanish. Using these 
matrix elements in Fermi's golden rule, Eq. (IT4t . gives a 
prediction of the tunneling rates 



FIG. 16. (Color online) Tunneling rates of regular states with 
quantum numbers n < 3 vs fc for a — 0.5 comparing the 
prediction of Eq. (|83p (solid lines) and numerical rates (dots) . 
The insets show the regular eigenfunctions iprcg^ and iprcg^ (as 
indicated by labels). In addition the asymptotic prediction of 
Eq. ((88)) is presented (dashed lines). 



mn / J 



X 3 



(82) 

The prime at the summation indicates that the sum over 
s excludes all multiples of three. The remaining integral 
can be solved analytically 86, Eq. 11.3.40], leading to 
the final result 



8 Jm+^-{jynnaf 



(83) 



This gives a prediction of direct regular-to-chaotic tun- 
neling rates of any regular state V'reg ^'^ chaotic sea in 
the mushroom billiard. The sum has its dominant con- 



tribution for s = 1 and evaluating Eq. ([S^]) up to s < 2 
gives sufficiently accurate predictions. 

It is worth to remark that a very plausible estimate 
of the tunneling rate is given by the averaged square of 
the regular wave function on a circle with radius a, i.e. 
the boundary to the fully chaotic phase space, yielding 
7mn = Nl^nJyniiranaY /2. Surprisingly, it is just about 
a factor of 2 larger for the parameters we studied. In 
Ref. 31] a related quantity is proposed, given by the 
integral of the squared regular wave function over the 
quarter circle with radius a. This quantity, however, is 
too small by a factor of order 100 for the parameters 
under consideration. 



2. Companson with numerical rates 

The eigenvalues and eigenfunctions of the mush- 
room billiard are determined by numerically solving the 
Schrodinger e quat ion. The improved method of partic- 
ular solutions |3ll . W\ allows a determination of the en- 
ergies E with a relative error sa lO"^"*. We analyze the 



widths Ai^ch.mn of avolded crossings between a given reg- 
ular state and typically 30 chaotic states under variation 
of the height I of the stem, starting with I — 0.3. From 
Eq. ([52)1 . 7 w (]A£'ch,mnl^ Ah)/8, we deduce the tunnel- 
ing rate where we use Acii — la- 

m 



[arcsin(a) -I- a^l — a?] /2 
as derived in Appendix [C] Note, that some pairs 
of regular states are very close in energy, e.g., i?20,i — 
£■16.2 ~ 10~^, such that their avoided crossings with a 
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chaotic state overlap, making a numerical determination 
of the smaller tunneling rate unfeasible within the pre- 
sented approach. 

Figure [16] shows the numerical tunneling rates 7mn for 
fixed radial quantum number ii = 1,2,3 and increasing 
azimuthal quantum number m for a — 0.5. It is com- 
pared to the theoretical prediction, Eq. which is 
connected for fixed n and increasing m by straight lines, 
giving rise to an apparently smooth curve. We find ex- 
cellent agreement for tunneling rates 7m„ over 18 orders 
of magnitude. The small oscillations which appear in the 
numerical rates on top of the exponential decay might be 
related to the two-level approximation for avoided cross- 
ings which we use for the numerical determination of the 
tunneling rates. 

To further test the prediction we determine the tun- 
neling rate of the regular state (m, n) = (30, 1) under 
variation of the stem width a. The results presented in 
Fig. [17] show a decrease of this tunneling rate which ap- 
pears faster than exponential with 1 — a. Again we find 
excellent agreement to numerical rates. Note, that the 
accuracy of the numerical method used for determining 
eigenenergies of the mushroom is best for a w 0.5 and 
declines for larger or smaller a. 

Another interesting question is how the tunneling rates 
from a given classical torus behave. For this we con- 
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FIG. 17. (Color online) Tunneling rates of the regular state 
ipreg^ vs the stem width a. We compare the prediction of 
Eq. (|83|l (solid lines) and numerical rates (dots). The insets 
show the regular eigenfunction V'reg^ at a = 0.34 and a = 0.74. 



FIG. 19. (Color online) Tunneling rates of regular states with 
quantum number m = 70 vs for a = 0.5 comparing the 
prediction of Eq. (|83p (solid lines) and numerical rates (dots). 
The insets show the regular eigenfunctions V'reg^ and V'reg^^- 



3. Approximation 



sider a sequence of regular states (m, n) which semiclas- 
sically locaHze on a torus characterized by an angular 
momentum pt > a. For each n we choose m such that 
Pmn — 'm/jnm ~ Pt- The resulting behavior of the tun- 
neling rates is presented in Fig. [TS] for pt = 0.6 and 
Pt = 0.8. A comparison of these predictions to numerical 
rates shows excellent agreement. 

For fixed azimuthal quantum number m and increasing 
radial quantum number n the tunneling rates increase. 
This behavior is presented in Fig. 1191 Again we find 
good agreement between the predictions of Eq. (|83l) and 
numerical rates. 
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FIG. 18. (Color online) Tunneling rates of regular states local- 
ized closest to a classical torus of angular momentum pt — 0.8 
and Pt — 0.6 vs k for a = 0.5. We compare the prediction of 
Eq. (|83p (solid lines) and numerical rates (dots). The insets 
show the regular eigenfunctions Vreg^ close to pt = 0.8 and 
ip^og'^, Vrcg'^5 close to Pt = 0.6. 



Let us now approximate Eq. for n and large 
wave numbers k in order to understand the exponen- 
tial behavior of the tunneling rates, which is visible in 
Figs. [TBI and [TBI with increasing k. First we consider the 
numerator of the leading term s = 1 in Eq. (|83p and 
use Ref. [S^, Eq. 9.1.63] for non-integer arguments of the 
Bessel function 



■{jmnO) < 



cxp I I m 



|)v/T-a: 



2 

'mn 



(l + Vl - "'ran) 



(84) 

with Qmn = jmnO'/i'm + 2/3). Equation (l84l) provides an 
upper bound of Jm+2/z{jmnO')- Numerically it has been 
confirmed, that a good approximation is given by this 
bound divided by to"^/^. 



1 



1 / ^ exp {mhmn) \ 



— 2" exp [ m 

m3 



bmn - 111 



1 + 6„ 



(86) 



where to = m -1- 2/3 and bmn = ■\/l — imn- The denomi- 
nator of the term s = 1 in Eq. (|83p can be approximated 
for n = 1 using Ref. 86, Eq. 9.5.18] 



Jni — li^jmn) — JfYii^jrnn) 



-1.1131TO-3 



TO3 



(87) 



We thus obtain for the tunneling rates assuming that 
Eq. ((87|) also approximately holds for n > 1 



7mn « - exp 2to 

TT 



In 



1 -|- hjji'ii 



(88) 



For fixed radial quantum number n and increasing az- 
imuthal quantum number m the tunneling rates decay 
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exponentially with k (X m. Figure [16] shows the com- 
parison to Eq. (dashed lines). We find agreement 
with deviations smaller than a factor of two. The pre- 
diction, Eq. (1551) . has a similar form as Eq. ([5^ which 
has been obtained for a quantum map with a harmonic 
oscillator-like regular island. This reflects the similarity 
of this system to the mushroom billiard. 



B. Annular billiard 

We consider the desymmetrized annular billiard char- 
acterized by the radius i? = 1 of the large semi-circle, 
the radius a of the small semi-circle, and the displace- 
ment w of this semicircle, see Fig. [^Ofa). In contrast to 
the sharply separated regular and chaotic dynamics in 
the mushroom billiard, the phase space of the annular 
billiard is more subtle. Again we find regions of reg- 
ular whispering-gallery motion. In the chaotic region, 
however, additional regular islands and partial barriers 
may be located, depending on the choice of a and w. 
This structure leads to the existence of so-called beach 
states which resemble regular states but are localized in 
the chaotic sea close to the border of the regular region. 
These beach states were described by Doron and Frischat 
[T^llTj who also studied the dynamical tunneling process 
in annular billiards. Their prediction of tunneling rates 
required fitting with a free parameter. In this section we 
want to apply the fictitious integrable system approach 
to find a prediction of direct regular-to-chaotic tunnel- 
ing rates which describe the decay of whispering-gallery 
modes into the chaotic sea. 
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FIG. 20. (Color online) (a) Schematic picture of the desym- 
metrized annular billiard showing the two coordinate systems 
used in the theoretical derivation of the direct regular-to- 
chaotic tunneling rates, (b) Poincare section at the semi-circle 
boundary (relative tangential momentum p vs arclength s) 
showing regular and chaotic regions for R = 1, a — 0.15, and 
w — 0.45 with illustrations of trajectories. 



In order to determine these direct regular-to-chaotic 
tunneling rates we proceed similar to the case of the 
mushroom billiard. We have to evaluate Eq. (fT3|) for 
the coupling matrix elements Vch,mn and then use Fermi's 
golden rule, Eq. (IT4l) . to determine the tunneling rates. In 
the first step the fictitious integrable system i/rcg and its 
eigenstates ip^g have to be defined. A natural choice for 
-fffog is the semi-circle billiard which exactly reproduces 
the whispering-gallery motion in the annular billiard. Its 
eigenstates are given by 

ip^g{r, ip) = 7V„„ J„ (j™„r) sin(m(/?) (89) 

in polar coordinates (r, (ys), where Jm denotes the mth 
Bessel function, jmn is the nth root of Jm, and Nmn = 
yjA/ir/ J,n-i{jmn)- The regular states are character- 
ized by the radial quantum number n = 1,2,... and 
the azimuthal quantum number m ~ 1,2,... . Hence, 
the tunneling rates describing the decay of the regular 
state V'reg(^i¥') ■^ill be labeled by jmn- Note, that for 
the annular billiard only those regular states semiclassi- 
cally exist which localize on tori with angular momentum 

Pmn = m/jmn > W + a. 

Evaluating Eq. (|13p in order to determine the cou- 
pling matrix elements Uch,mn between the regular and 
the chaotic states, an infinite potential difference arises 
within the small disk of radius a between the Hamilto- 
nian H of the annular and i?rcg of the semi-circle billiard, 
H — i/reg = oo. At the same time for the chaotic states 
"^ch — holds in that region, which leads to an unde- 
fined product "oo • 0" . Similar to the approach presented 
for the mushroom billiard we circumvent this problem 
by considering a finite potential difference W for which 
at the end the limit — >■ oo is performed. In contrast 
to the mushroom billiard the area of the annular billiard 
n is included in the area of the semi-circle billiard iljcg, 
O C rirog- We find that in this case the derivative of the 
chaotic states ipch{q) enters in Eq. (fT3t . We obtain 

VcKrnn = a j Vj™g"(a, d) dp^P = «, ^) (90) 


= aNmn J d-dJmijmnr) sm{m(p)dpipchia, ■&) , (91) 



where we introduce polar coordinates (p, d) such that 
X = pcos{'d) + w and y = psin(i9). Only the integration 
over 1? from = to = tt along the small semi-circle 
of radius p — a remains. Note, that the regular eigen- 
functions are given in polar coordinates {r{p, i?), (p{p, "&)), 
while we integrate along p = a, see Fig. BUT a'). Along this 
half-circle the regular wave function is largest near the 
point (p, 1?) = (a, 0) for w > 0. Hence, a random wave 
description for the chaotic states ipchip, t^) bas to respect 
the Dirichlet boundary conditions on the line y = and 
on the small semi-circle of radius a. 

Such a random wave model can be constructed using 
the solutions of the annular concentric billiard [88), §25] 
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as base functions in polar coordinates (p, i?) 



: Cs sin(s'!9) 



s=l 



J,{kp)Ys{ka) ~ Js{ka)Y,{kp) 
^Js{kaY+Ys{kaY 



(92) 



in which the Cg are Gaussian random variables with 
(cs) = and (cgCt) — Sgt- Furthermore, Yg denotes the 
sth Bessel function of the second kind. The normaliza- 
tion constant 2 has been obtained numerically, such that 
(IV'chP) = l/-4ch holds far away from the boundary. The 
chaotic states defined by Eq. ((92|) fulfill the Schrodinger 
equation at arbitrary energy E — k^ . Similar to the 
mushroom billiard we do not require that the chaotic 
states decay into the regular island, as Eq. (|9T|) is an 
integral along a line of the billiard which is not hit by 
any regular whispering-gallery trajectory. Near the hor- 
izontal boundary and away from the small circle (iV'chP) 
recovers the behavior 1 — Jo(2fc|p— a\) 84, 85]. Using 
the radial derivative of i}}c\i{p, "&) at p = a one finds 



dpipchip = 0,1?) 



2k 

Cs sm{s'd)Rs{ka), 



(93) 



in which we use Js+i{ka)Ys{ka) — Js{ka)Ys+i{ka) 
2/{-Kka) m, Eq. 9.1.16] and introduce 



Rs{ka) 



nkay/ Js{ka)^ + Ys{kay 



(94) 



With this result we calculate the tunneling rates at en- 
ergy E,nn = jfnn 

oo 



the regular whispering-gallery modes remain almost un- 
changed. Note, that under this parameter variation the 
phase-space structure in the chaotic region changes from 
macroscopically chaotic to a situation in which additional 
regular islands and partial barriers appear. With increas- 
ing k, however, a smaller parameter variation is required 
and these problems become less relevant. 

Figure [^l] shows the numerical tunneling rates for 
w + a — 0.6 and the analytical prediction, Eq. (^5]) . 
for a — 0.15 and w — 0.45. Variations of a with con- 
stant w -|- a, as necessary for the numerical rates, affect 
the prediction by at most a factor of two (not shown). 
The comparison in Fig. [5T] gives qualitative agreement. 
While some tunneling rates agree with the prediction, 
deviations of a factor of 100 appear for other rates. 

We believe that these deviations are artifacts of our 
numerical procedure to determine tunneling rates from 
avoided crossings They occur most likely due to 

beach states which exist in the chaotic region of phase 
space and look similar to regular states though no quan- 
tizing tori can be associated with them [ll,[i3. Numer- 
ically we need to analyze avoided crossings between the 
regular mode and all modes outside the regular island 
under variation of the shape of the billiard. However, as 
beach states almost behave like regular states, their en- 
ergy only slightly varies if the boundary of the billiard is 
changed. Hence, they rarely show avoided crossings with 
the regular state. In addition the chaotic states concen- 
trate further away from the regular island and thus give 
rise to considerably smaller avoided crossings. This leads 
to artificially reduced numerical tunneling rates. It would 
be desirable to determine numerical tunneling rates for 
the annular billiard by other means, e.g., by opening the 
billiard, to obtain a more quantitative verification of the 
prediction Eq. ([55)) . 



M sin(sz9) Jm(jm„r(i9)) sm{mip{d)) 



(95) 



This is our final result predicting the decay of a regular 
state ■(/'™g located in the regular whispering-gallery re- 
gion to the chaotic sea. In Eq. (^5]) . in contrast to the 
result for the mushroom billiard, Eq. (l83l) , the term s = 1 
is typically not the most important contribution. Here 
the dominant s increases with energy. In contrast to the 
prediction of Refs. [l^ [i3 no fitting is required. 

The numerical determination of tunneling rates using 
avoided crossings is more difficult for the annular bil- 
liard than in the case of the mushroom billiard: In order 
to affect the chaotic but not the regular component of 
phase space under parameter variation we increase the 
radius a of the small inner circle and move its center 
position w such that its rightmost edge at w -I- a is con- 
stant. We increase a starting from a = 0.05 until 30 
avoided crossings have occurred or a = 0.4 is reached. 
This procedure drastically affects the chaotic states while 
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FIG. 21. (Color online) Tunneling rates from regular states 
with quantum numbers n < 3 vs fc for w + a — 0.6 comparing 
the prediction of Eq. (|95p (solid lines) and numerical rates 
(dots). The insets show the regular eigenfunctions 'i/'res^ a^nd 

reg ■ 
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C. Disordered wire in a magnetic field 



We study two-dimensional nanowires [55|, [561, l89[ with 
one-sided disorder, see Fig.[22la). in the presence of a ho- 
mogeneous magnetic field B perpendicular to the wire. 
Such nanowires have a mixed phase space, see Fig.[22][b). 
Orbits which only hit the lower flat boundary are regu- 
lar skipping orbits while those which are reflected at the 
upper disordered boundary numerically show chaotic mo- 
tion. The phase space of the wire has a sharp transition 
from regular to chaotic dynamics. There are no resonance 
chains inside the regular island and there is no relevant 
hierarchical region. 

In Refs. [55, . .56] it was shown that in such wires the 
localization lengths ^ increase exponentially with increas- 
ing Fermi wave number kp of the electrons. This can be 
explained by the mixed phase-space structure giving rise 
to dynamical tunneling between the regular island and 
the chaotic sea. The dynamical tunneling rates 7 are di- 
rectly related to the localization lengths, ^ cx I/7. We 
apply the fictitious integrable system approach to deter- 
mine localization lengths ^ and compare the results to 
the analytical prediction derived in Refs. [ssHSGf. 

The nanowire with a disordered boundary is assembled 
from L rectangular elements to which two leads of width 
W are attached. The nth element has a length / — W/5 
and a height h„ which is uniformly distributed in the in- 
terval [T^-^/2, W+S/2] with 5 = 2W/3. The heights /i„ 
are allowed to take M — 20 different values including the 
boundaries of the interval. For increasing values of the 
Fermi wave number kp the magnetic field B is adjusted 
such that the cyclotron radius Vc — KkpjieB) remains 
constant, using = 3W. This leaves the classical dy- 
namics unchanged and the semiclassical limit is given by 
kp — > 00. 
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FIG. 22. (Color online) (a) Schematic picture of a mag- 
netic wire with one-sided disorder to which leads of width 
W are attached. The wire is composed of rectangular el- 
ements of length I and a uniformly distributed height in 
\W — & /2,W -\- 5 /2\. For the tunneling process contributions 
(i) along the vertical and (ii) along the horizontal parts of the 
upper boundary are relevant, (b) A Poincare section {y,Py) 
at fixed X shows a large regular island (red lines) and the 
chaotic sea (blue dots). 



In order to apply the fictitious integrable system ap- 
proach we have to consider a closed billiard. It is ob- 
tained by imposing periodic boundary conditions at the 
positions where the leads are attached. The tunneling 
rate is given by Fermi's golden rule, Eq. (fT4)) . using the 
coupling-matrix elements between regular and chaotic 
states, Eq. (fTS]) . Similar to the examples discussed previ- 
ously we have to find a fictitious integrable system -fficg 
which resembles the regular part of the phase space and 
extends it beyond. A reasonable choice is a billiard of 
length LI with periodic boundary conditions at a; = 
and with a boundary at y = but open for y > 0. In 
a magnetic field it shows exactly the same regular phase 
space as the wire with one disordered boundary and is 
completely integrable. Its eigenfunctions tpl!^g{x,y) are 
given by [90| 



1 



Ak-mX 



Zm{y), 



(96) 



where Zm{y) solves a one-dimensional Schrodinger equa- 
tion with an effective potential due to the magnetic field 
and m > 1 is the quantum number in transversal di- 
rection. Using this regular system and its eigenstates in 
Eq. ([Ta]) gives 



ds dNA^.{x{s),v{s)) iJ;::Jx{s),yis)). (97) 



Vch,; 



As for the annular billiard the derivative of the chaotic 
states dNipch{x{s),y{s)) is in the normal direction along 
the disordered boundary C, which is parameterized by s. 

The most important contributions of the integral in 
Eq. (P?]) are from rectangular elements of the wire which 
have the smallest height, i.e. /i„ = W — S/2. This follows 
from the exponential decay of the regular wave functions 
ip™g^ in y-direction. There are L/AI such rectangular el- 
ements. The upper boundary of the wire near such an 
element is composed of (i) two vertical parts (x = nl, y € 
[h„-i,W-S/2] and x = {n+l)l, y e [W-S/2, hn+i]) and 
(ii) one horizontal part {x S [nl, (n + 1)Z], y — W — S/2). 

(i) Let us first consider one vertical contribution to the 
coupling-matrix elements Vch,m for the fixed ^-coordinate 
X* = nl. For the chaotic states tpch we employ a random 
wave model with wave length kp such that the Dirichlet 
boundary condition at x = x* is fulfilled 



:^^CsJs{kpr) sin 



Sip - 



,(98) 



where the polar coordinates (r, ip) are defined hy x = 
X* + r cos (p and y — W — S/2 + r sin tp. The coefficients 
Cs are Gaussian random variables with mean zero and 
(csCj) = 6 St- It will turn out that the form of Eq. (|98l) is 
more convenient for the following evaluation than a plane 
wave ansatz. The derivative of '0ch with respect to x at 
X = X* gives 



dxi^chix = x*,y) 



^,/s{kpy)s^ (99) 
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where y = y — {W — 5/2). Hence, the contribution of 
one vertical part of the boundary to the coupling-matrix 
element is (up to a phase e''^™'^ ) 



Vch.'. 



Js{kFy)s 




Z.^{y) dy.(lOO) 



Here we have replaced the upper integration limit by 
oo as Zm{y) decays exponentially. The sum of 2L/M 
such coupling matrix elements gives according to Fermi's 
golden rule, Eq. ([T^ . the direct regular-to-chaotic tun- 
neling contribution of the vertical boundaries 



v(0 



Ml 



Jsikpy) 



V 



Z^{y)dy\ , (101) 



where we assume independent coefficients Cg for each ver- 
tical boundary. 

(ii) For the horizontal boundaries we consider the regu- 
lar wave functions ipl'^g{x,y = 0) = ^o^Qik^{ni+x} ^ where 

X = X — nl and Z^^ — Zm{y — 0). For the chaotic wave 
function a random wave model respecting the Dirichlet 
boundary at y = is used. 



V'ch(i,y) = 



NA, 



ch 



N 

E 

s=l 



{{kpx cos i?s + V^s ) 



ce'^'"'' " , ^s, sin(fci?y sini^s), 

(102) 

where the coefhcients Cg are Gaussian random variables 
with mean zero as well as {cgCt) — 5st and the angles "ds 
and are uniformly distributed in [0, 27r). The deriva- 
tive of '0ch with respect to y at y = reads 



dyii)ch{x,y = 0) = 



ch 



N 

E 

s=l 



csini^, gi(fcf£cosi;,+v,)^ 

Hence, we obtain for the contribution of one horizontal 
part of the boundary to the coupling-matrix element (up 
to a phase e'*''™"') 



^ch.m 1—1 T 1 / 







(104) 

The sum of L/M such coupling matrix elements gives 
according to Fermi's golden rule, Eq. p^ . the direct 
regular-to-chaotic tunneling contribution of the horizon- 
tal boundaries 



2nMl 



M sin^ ■& 



2-2cos{lkp cos-d+^ 



(cos ^+17) 



(105) 

where we assume independent coefficients Cg for each hor- 
izontal boundary. The total tunneling rate is given as the 
sum of the two contributions, Eq. (|10ip from the vertical 



parts of the boundary and Eq. (llOSp from the horizontal 
parts of the boundary, 



7r- 



(106) 



Here we use the approximation that random wave models 
for the horizontal and vertical boundaries are indepen- 
dent. 

We now compare this result for the dynamical tunnel- 
ing rates 7„i to the semiclassical prediction of localization 
lengths for the wire with one-sided disorder derived in 
Ref. [S^ . For the connection between tunneling rates and 
localization lengths we use their inverse proportionality 
[9l|,[92. One finds 

where Vx,m is the velocity of the mth mode in x-direction 
and the localization length is measured in units of the 
length I of the rectangular elements. Figure [23l compares 
the localization lengths obtained by the fictitious inte- 
grable system approach using Eqs. (|106p and (|107p with 
numerically determined Zm (y) (solid lines) to the analyt- 
ical prediction of Ref. [5g (dashed lines) showing good 
agreement with deviations smaller than a factor of two. 

Note, that while in Ref. [H^ the localization length 
is Z-independent we find for our horizontal contribution 
an Z-dependence due to the integral in Eq. (|104p . This 
becomes relevant for I 3> W and can be corrected by 
either changing the upper limit of integration in Eq. (jl04p 
from I to 1/kp or by choosing an improved i/rcg- This 
accounts for the fact that far away from the corners of 
the disordered boundary there is no tunneling. 



D. Optical microcavities 

Optical microcavities in which photons can be con- 
fined in three spatial dimensions are a subject of inten- 
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FIG. 23. (Color online) Localization lengths of transversal 
modes in a magnetic wire with one-sided disorder for the mode 
numbers m < 4 vs krW/TT for l/W = 1/5 and S = 2W/3. We 
compare the prediction of Eqs. (I106P and (|107|) (solid lines) 
and the analytical result of Ref. 56] (dashed lines). 
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sive research as they are relevant for apphcations such as 
ultralow-threshold lasers. Especially whispering-gallery 
cavities such as microdisks, microspheres, and micro- 
toroids have been investigated as they can trap photons 
for a long time near the boundary by total internal reflec- 
tion at the optically thinner medium. The corresponding 
whispering-gallery modes have a very high quality factor 
Q. While the microdisk emits the photons isotropically, 
cavities with deformed surfaces may additionally lead to 
directed emission [13, [zl, IH, HI] . 

The quality factors Q of optical microcavities can be 
connected to dynamical tunneling rates 7. These rates 
can be determined using the fictitious integrable system 
approach [sj], which will be summarized in this section. 
For a mode in an open cavity the quality factor Q is re- 
lated to the corresponding resonance with complex wave 
number k = Re(fc) + ilm(fc) via 



Q 



Re(/c) 
' 2Im(fc) ' 



(108) 



For cavities with a mixed phase space the quality factor 
Q of a regular mode has two contributions 



1 



1 



'.Id 



1 



Qdyn 



(109) 



Here Qid accounts for the coupling of the regular mode 
to the continuum. In the case of the circular microcavity 
it is described by a quasi one-dimensional barrier- 
tunneling process through an angular momentum barrier. 
It can be predicted by means of WKB theory or using 
the numerically determined decay rates of the modes in 
the circular cavity. Note, that for this contribution the 
mixed phase-space structure is irrelevant. The second 
contribution, Qdyn, is given by dynamical tunneling from 
the regular mode to the chaotic sea, which is strongly 
coupled to the continuum. Here we assume that there are 
no further phase-space structures within the chaotic sea 
that affect the quality factor. A priori it is not obvious, 
which of the two contributions will dominate. 

In order to determine Qdyn the fictitious integrable sys- 
tem approach can be employed in analogy to hard-wall 
billiards. It is extended to open cavities in the follow- 
ing way: (i) As a fictitious integrable system -ffrcg one 
chooses a cavity such that it resembles the regular dy- 
namics of H. The quantum system has resonance states 
^reg- (ii) As a model for the chaotic resonances il^ch a- 
random wave model is used, which in addition fulfills the 
relevant cavity boundary conditions, (iii) The tunneling 
rate 7 determines the quality factor Qdyn by 



2Rc(fc)2n2 



7 



(110) 



where no is the refractive index of the cavity. Here it 
is used that the imaginary part of the wave number is 
connected to the decay rate of the resonant state via 
7 = -4n^Re(fc)Im(fc) and Eq. (ITU5)) . 



This approach was developed and applied to the annu- 
lar microcavity, a microdisk with an air hole, in Ref. j34| . 
It has a particularly interesting geometry which allows for 
unidirectional emission and high quality factors simulta- 
neously. An analytical expression for the quality factor 
Qdyn was derived which is in very good agreement with 
the full numerical simulations of Maxwell's equations. 



VI. SUMMARY AND OUTLOOK 

We study the direct regular-to-chaotic tunneling pro- 
cess in systems with a mixed phase space. It is dominant 
in the regime, hcs ^ ^reg, where fine-scale structures of 
the phase space such as nonlinear resonances are not re- 
solved by quantum mechanics. To describe this tunneling 
process we introduce the fictitious integrable system ap- 
proach. It uses a decomposition of the Hilbert space into 
two parts which account for the regular and the chaotic 
dynamics. This leads to a formula which predicts the di- 

to 



rog/ 



rect tunneling rate 7,„ of the mth regular state 
the chaotic sea, Eq. ([7]). The fictitious integrable system 
has to be chosen such that its dynamics resembles the 
regular dynamics of the original mixed system as closely 
as possible and extends it beyond its regular region. The 
determination of the fictitious integrable system is the 
most difficult step in the application of this approach. 
For maps the Lie transformation or methods based on 
the frequency map analysis can be used. A general pro- 
cedure for billiards is under development. 

In Sec. IIVI dynamical tunneling rates are determined 
numerically and compared to the prediction of the ficti- 
tious integrable system approach for different quantum 
maps. We find excellent agreement over several orders 
of magnitude in 7. For a harmonic oscillator-like island 
embedded in a chaotic sea it is possible to derive a semi- 
classical expression for the tunneling rates which depends 
on the area of the regular region and the effective Planck 
constant only. Furthermore, in Sec. |V] we apply the ap- 
proach to billiard systems, where we use random wave 
models to describe the chaotic states. For the mushroom 
billiard the prediction is in excellent agreement with nu- 
merical and experimental data. Finally, we apply the 
approach for direct regular-to-chaotic tunneling rates to 
the annular billiard, nanowires with one-sided disorder 
in a magnetic field, and to the quality factors of optical 
microcavities. 

In the semiclassical regime, hcs <C ^rog, nonlin- 
ear resonances lead to enhanced tunneling rates due 
to resonance-assisted tunneling. This resonance-assisted 
tunneling mechanism is combined to direct regular-to- 
chaotic tunneling in Ref. [sH], leading to a prediction 
of tunneling rates which is valid from the quantum to 
the semiclassical regime. For this prediction the direct 
tunneling rates discussed in this paper, which can be de- 
termined using the fictitious integrable system approach, 
are essential. 

Ultimately, it is the aim to obtain a complete semiclas- 
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sical treatment of the approach also for systems with a 
generic regular island. This theory should predict tun- 
neling rates and their dependence on the effective Planck 
constant /log directly from classical properties of the reg- 
ular island, such as its size, shape, and winding num- 
ber and properties of the chaotic sea, such as unstable 
fixed points, partial barriers, and transport properties. 
A promising approach seems to be the complex-path for- 
malism which has been successfully applied to study the 
tunneling tails of a time-evolved wave packet in mixed 
regular-chaotic svstems |12h14|. Alternatively, a com- 
plex time approach [23.l25l| can be considered. Using such 
semiclassical approaches for the prediction of regular-to- 
chaotic tunneling rates should give further insight into 
the dynamical tunneling process. 
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Appendix A: Dimensionless Fermi's golden rule for 
time-periodic quantum systems 

In order to predict dynamical tunneling rates 7m we 
use Fermi's golden rule for systems with a discrete spec- 
trum. This rate describes the decay e^'*'™' of the mth 
regular state to the chaotic sea at most up to the Heisen- 
berg time th = /loff/Ach- Com par ed to the standard 
derivation of Fermi's golden rule [96| for the continuous 
case, the incoherent integral over the continuum states is 
replaced by the incoherent sum over the discrete chaotic 
states. One finds 



27r 



(IV^c 



ch. 



(Al) 



with the chaotic density of states pch and the average 
of the modulus squared of the coupling matrix elements 
Ki:h,m of the Hamiltoniau between energetically close-by 
chaotic states and the mth purely regular state. 

We now apply Eq. (jAll) to time-periodic systems with 
period r and quasi-energies in the interval [0, Hlo] with 
w = 27r/T such that pch = A"ch/(^)- We introduce di- 
mensionless quantities 7m = FmT and Wch,m = K;h,mT/ft 
leading to 



7m = Nch{\Vch,: 



E 

ch 



(A2) 



Here Vch,m = (V'chjt^lV'rcg) is the coupling matrix element 
of the time evolution operator U over one period and the 
mean level spacing is 2tt/N . 



Appendix B: Adapted eigenstates of the harmonic 
oscillator 



For the maps 2?, which for i? = show a tilted and 
squeezed harmonic oscillator-like regular island, purely 
regular states tp™g{q) can be constructed by tilting and 
squeezing the eigenfunctions of the harmonic oscillator 
accordingly. One finds [o^] 



1 /Re(cr) 



H„ 



'Re(cr) \ 
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with the Hermite polynomials H„i and the complex 
squeezing parameter ct, which can be obtained from the 
tilting angle 6 of the elliptic island with respect to the 
momentum axis and the ratio uj of its half axes with 



Re(cr) = - 
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icos2(6') -hwsin2(6l)' 
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lm{(7) ^ Re{(7) y ujj sm{0) co8{0). (B3) 

Here 9 and ui can be determined from the linearized dy- 
namics of the classical map around the fixed point. 
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where c = y/{Mi2 + M21Y + {M22 - MnY and 
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Appendix C: Derivation of for the mushroom 
billiard 



We want to show how the area v4ch, which is the area of 
the billiard times the fraction of the chaotic phase-space 
volume, is determined for the mushroom billiard. 

All trajectories which enter the stem of the mushroom 
or cross the half-circle r = a in the cap are chaotic irre- 
spective of their momentum direction. Additionally the 
region r > a contains some chaotic dynamics. We denote 
its contribution to Ach by / and obtain for the mushroom 
billiard 



A, 



cli 



2la- 



I. 



(CI) 
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FIG. 24. (Color online) All trajectories of the mushroom bil- 
liard passing through the point (r, if) with r > a (dot) and 
angle < i9c cross the half-circle of radius a (dashed line) 
and are chaotic. For angles |i? -I- 7r| < i^c the time-reversed 
trajectories cross the half-circle and they are chaotic as well. 



In order to determine / we consider the fraction P{r) of 
chaotic trajectories in the cap for r > a. It depends on 
the radial coordinate r only 



2 /a\ 
P{r) ~ — arcsin I — 



(C2) 



as can be seen from Fig. [24l 

Integrating over the region r > a in the cap of the 
mushroom gives 



IT R 



j J rdr dipP{r) 



(C3) 



arcsin f + a^/R^ - - iyra^. (C4) 
\R/ 2 



Inserting this expression into Eq. (|Cip we finally obtain 
^ch = 2^a 



R^ arcsin yj^j + ay/R^—a^ 



(C5) 



as the the area of the mushroom billiard times the frac- 
tion of the chaotic phase-space volume. For the desym- 
metrized mushroom used in Sec. IV Al this result has to 
be divided by two. 
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